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SUMMARY 


This  report  describes  studies  of  CW  laser-heated  rocket  propulsion, 
in  both  the  flowing-core  and  stationary-core  configurations,  with  most 
attention  focused  on  the  former.  (It  is  the  second  such  study  performed 
by  PSI  following  the  work  described  in  Ref.  1).  In  the  present  work,  the 
laser  radiation  considered  was  10.6  pm,  and  the  working  gas  was  unseeded 
hydrogen.  The  areas  investigated  included  initiation  of  a  hydrogen  plasma 
capable  of  absorbing  laser  radiation,  the  radiation  emission  properties  of 
hot,  ionized  hydrogen,  the  flow  of  hot  hydrogen  while  absorbing  and  radiating, 
the  heat  losses  from  the  gas  and  the  rocket  performance.  The  stationary- 
core  configuration  was  investigated  qualitatively  and  semi-quantitatively. 

It  was  found  that  the  flowing-core  rockets  can  have  specific  impulses 
(Igp)  between  1500  and  3300  sec.  They  are  small  devices,  whose  heating 
zone  is  only  a  millimeters  to  a  few  centimeters  long,  and  millimeters  to 
centimeters  in  radius,  for  laser  power  levels  varying  from  10  to  5000  kW, 
and  pressure  levels  of  3  to  10  atm.  Heat  protection  of  the  walls  is  a  vital 
necessity,  though  the  fraction  of  laser  power  lost  to  the  walls  can  be  as  low 
as  10%  for  larger  powers,  making  the  rockets  thermally  efficient. 

A  number  of  major  areas  of  uncertainty  have  been  identified  for 
further  exploration.  Chief  among  these  is  the  properties  of  laser  supported 
combustion  (LSC)  waves  in  hydrogen.  Others  are  an  efficient  method  of 
heat  protection,  the  effect  of  radial  temperature  profiles  on  radiation  loss, 
and  the  effect  of  wider  variation  of  operating  parameters.  For  the  stationary- 
core  concept,  absorption  lengths  of  hot  hydrogen  radiation  in  cold  hydrogen 
are  needed,  and  a  study  of  entrainment  of  the  core  gas  by  the  flowing 
propellant  gas  should  be  made. 
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1.  INTRODUCTION 


This  is  the  final  report  on  the  second  contract  issued  by  NASA/ 

Lewis  Research  Center  to  Physical  Sciences  Inc.  for  a  study  of  a  laser- 
heated  rocket.  The  results  of  the  first  contract  were  reported  in  Ref.  1. 

The  present  study  dealt  with  the  details  of  a  number  of  aspects  of  a 
laser-powered  rocket  engine  absorbing  CO2  laser  radiation  (10.6  p.m)  by 
inverse  Bremsst  rah  lung,  and  using  hydrogen  as  the  working  gas.  There 
were  two  general  areas  of  study.  The  first  was  the  initiation  of  a  high 
temperature  plasma  capable  of  absorbing  the  laser  radiation.  The  second 
was  the  analysis  of  a  steady  coaxial-flow,  laser-powered  rocket  engine  in 
two  different  configurations.  One  configuration  was  a  flowing  core  of  plasma, 
absorbing  the  laser  energy,  being  heated,  and  expanding  to  generate  thrust. 
This  core  is  surrounded  by  a  cooling  buffer  gas  flowing  the  annulus  between 
the  core  and  the  engine  walls.  The  other  configuration  was  a  stationary 
stable  plasma  absorbing  the  laser  energy,  and  transmitting  it  to  a  propel¬ 
lant  gas  flowing  around  it,  which  in  turn  expands  to  generate  thrust. 

The  general  geometries  of  these  two  configurations  are  shown  in 
Fig.  1.  1.  In  the  flowing-core  concept,  cold  hydrogen  flows  into  the  engine, 
parallel  to  the  laser  beam,  and  flows  into  a  standing  laser  supported 
combustion  (LSC)  wave.  This  wave  heats  the  hydrogen  up  to  a  temperature 
at  which  it  can  absorb  laser  energy  by  inverse  Bremsstrahlung.  The  hydro¬ 
gen  is  further  heated  by  this  process  while  flowing  in  a  variable  area 
streamtube,  until  all  the  laser  energy  is  absorbed.  The  temperature  of  the 
hydrogen  is  sufficiently  high  that  it  is  a  dissociated  and  ionized  plasma  behind 
the  LSC  wave,  and  radiates  energy.  When  all  the  laser  energy  is  absorbed, 
the  heated  plasma  expands  through  a  nozzle  to  generate  thrust.  This 
plasma  core  streamtube  will  transmit  energy  through  its  side  walls  by  both 
convection  and  radiation,  and  to  protect  the  rocket  walls  a  co-flowing  buffer 
gas  stream  not  heated  by  the  laser  may  be  necessary.  The  amount  of  buffer 
gas  flow  necessary  is  determined  by  the  amount  of  energy  lost  radially  by 
the  plasma  core. 

In  the  stationary- core  concept,  a  hot  core  of  stationary  hydrogen 
continually  absorbs  laser  radiation  and  re-radiates  it  to  a  surrounding 
annular  streamtube  of  flowing  hydrogen  propellant  gas.  The  propellant  gas 
absorbs  the  core  radiation  and  becomes  heated.  After  passing  the  core,  it 
is  expanded  to  generate  thrust.  The  propellant  gas  also  acts  as  a  buffer  to 
protect  the  walls  from  the  hot  core. 
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Clearly,  these  kinds  of  engines  depend  on  a  large  number  of  physical 
processes.  There  is  laser  radiation  absorption,  for  which  the  inverse 
Bremsstrahlung  mechanism  has  already  been  defined  as  the  desirable 
mechanism  in  Ref.  1.  There  is  the  radiation  emission  by  the  hydrogen 
plasma,  which  is  a  loss  in  the  flowing-core  concept  and  the  main  energy 
transfer  mechanism  in  the  stationary- core  concept.  There  is  the  absorption 
of  plasma  radiation  by  colder  hydrogen,  which  heats  the  propellant  gas  in 
the  stationary- CO  re  concept,  and  the  on-coming  gas  in  the  LSC  wave.  This 
absorption  also  can  protect  the  walls.  There  is  convective  heat  transfer 
from  the  hydrogen  plasma  core  to  the  surroxinding  environment.  There  is 
the  flow  of  heated  hydrogen  in  a  variable  area  streamtube  while  dissociating 
and  ionizing.  There  is  the  LSC  wave  mechanism  by  which  a  laser  can 
convert  room-temperature  hydrogen  to  a  plasma  which  can  absorb  the  laser 
energy.  There  is  the  interfacial  mixing  as  the  propellant  gas  flows  past  the 
core  in  the  stationary- core  concept. 

In  addition  to  these  processes  important  to  the  steady  state  operation 
of  the  rocket,  other  physical  phenomena  are  pertinent  to  the  initiation  of  the 
absorbing  plasma.  These  include  the  possible  mechanisms  of  gas  breakdown 
by  electric  discharge,  by  a  focussed  laser,  or  by  absorption  of  laser  energy 
by  molecules  or  particles. 

Although  all  these  processes  are  involved  in  the  operation  of  a  laser- 
heated  rocket,  not  all  of  them  are  studied  in  detail  in  this  report.  Here,  the 
emphasis  is  on  the  flowing-core  concept,  and  the  processes  important  to 
its  operation  are  the  ones  to  which  most  attention  is  given.  The  work  can 
be  viewed  as  an  extension  of  the  studies  presented  in  Ref.  1. 

Information  is  presented  on  the  three  mechanisms  of  plasma  initia¬ 
tion  by  electric  arc  discharge,  by  laser-induced  breakdown,  and  by  molecu¬ 
lar  and  particulate  abosrption,  in  Secs.  2,  3  and  4.  A  detailed  study  of  the 
radiative  emission  from  hydrogen  is  given  in  Sec.  5,  leading  to  expressions 
useful  for  coupling  the  emission  to  the  fluid  mechanics  of  flowing  hydrogen. 
Section  6  presents  an  estimate  of  the  intensity  vs  wave  speed  for  LSC  waves 
in  hydrogen  by  a  rough  extension  of  one  of  the  theories  for  air.  In  Sec.  7 
is  a  detiled  description  of  a  hydrogen  plasma  in  chemical  equilibrium, 
absorbing  laser  radiation  by  inverse  Bremsstrahlung,  radiating,  and  flow¬ 
ing  through  a  variable  area  channel.  Section  8  presents  calculations  on  the 
heat  loss  by  both  radiation  and  convection  of  the  flowing  plasma  core,  and 
an  estimate  of  the  radiation  which  can  be  absorbed  by  particles  in  a  buffer 
gas.  In  Sec.  9,  the  performance  of  the  flowing  plasma  core  as  a  rocket  is 
given,  based  on  the  description  of  the  core  in  Sec.  7.  The  stationary- core 
rocket  is  discussed  in  Sec.  10  in  a  qualitative  and  semi -quantitative  way, 
based  on  estimates  of  the  physical  phenomena  involved.  The  conclusions 
obtained  so  far  in  the  study  of  laser-heated  rockets  are  presented  in  Sec.  11. 
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The  cases  calculated  for  the  flowing-core  streamtube  were  laser 
powers  of  10,  100  and  5000  kW,  each  at  3  and  10  atm  initial  pressure. 

These  choices  were  made  to  provide  information  on  sizes  of  engines  in  which 
NASA/Lewis  has  expressed  interest.  An  experimental  engine  designed  to 
operate  at  10  kW,  3  atm  is  being  designed  and  built  by  Rocketdyne  Division 
of  Rockwell  International,  under  contract  to  NASA/Lewis,  and  the  same 
organization  is  doing  a  study  of  a  5000  kW  engine.  PSI  has  aided  Rocketdyne 
with  the  fluid  mechanics  and  radiation  aspects  of  their  program,  under  a 
subcontract.  The  choice  of  a  100  kW  laser  was  made  to  provide  an  inter¬ 
mediate  size  between  10  and  5000  kW,  which  might  some  day  be  tested.  The 
calculations  at  the  higher  pressure  of  10  atm  were  made  to  show  the  effect 
of  pressure  level  on  the  core  flow. 

In  addition  to  the  laser  power,  and  pressure  level,  there  are  a 
number  of  other  parameters  which  could  be  varied  in  the  calculations. 

One  is  the  laser  intensity  which  determines  the  mass  flow  rate  per  unit 
area  into  the  core.  In  the  present  work  this  was  held  fixed  at  3.  67E  9  W/m^ 
(3.  67E  5  W/cm^),  which  seems  a  value  at  which  the  ESC  wave  could  be 
supported.  When  better  knowledge  of  ESC  waves  in  hydrogen  is  available, 
other  values  may  be  of  interest.  (The  cross-sectional  area  of  the  rocket 
is  inversely  proportional  to  this  intensity  for  a  fixed  laser  power. )  A  ' 
second  parameter  which  could  be  varied  is  the  streamtube  velocity  distribu¬ 
tion,  which  was  also  held  fixed  in  the  present  calculations.  Such  parametric 
calculations  could  be  performed  by  a  slight  improvement  of  the  core  flow 
program  used  here,  and  would  aide  in  mapping  out  the  regions  of  interest 
for  laser-heated  rockets. 


-4- 


2.  PLASMA  INITIATION  BY  ELECTRIC  DISCHARGE 


One  possible  technique  for  initiating  laser  supported  combustion 
(LSC)  waves  in  hydrogen  gas  is  through  use  of  an  electric  discharge.  The 
purpose  of  the  discharge  is  to  supply  an  initial  density  of  electrons  in  the 
gas  so  that  the  degree  of  absorption  of  laser  energy  by  inverse  Bremsstrah- 
lung  is  sufficient  to  drive  the  gas  to  the  required  temperature  for  formation 
of  a  stable  LSC  wave.  Electric  discharges  have  previously  been  success¬ 
fully  used  to  initiated  LSC  waves  in  air^.  No  such  studies  are  available 
for  hydrogen.  Thus  the  first  question  to  be  answered  is  what  type  of 
discharge  is  best  suited  for  LSC  ignition. 

Electrical  discharges  can  be  broken  into  two  basic  groups,  depending 
on  whether  they  are  or  are  not  self-sustaining.  In  the  latter  group,  desig¬ 
nated  as  Townsend  discharges,  electrons  are  formed  at  the  cathode  by  a 
technique  such  as  thermionic  emission;  these  are  then  accelerated  by  the 
electric  field  and,  depending  upon  the  field  strength,  these  initial  electrons 
can  further  ionize  the  gas  within  the  discharge  gap.  The  characteristic 
electron  densities  in  such  discharges  are  ~  lO^Vcc  and  the  electron  temp¬ 
eratures  can  be  several  electron  volts,  well  out  of  equilibrium  with  the 
gas  translational  temperature. 

As  the  voltage  is  increased  across  the  electrodes  a  point  is  reached 
where  a  spark  is  created  across  the  gap.  This  spark  is  a  result  of  gas 
breakdown  and  will  occur  even  if  the  cathode  does  not  emit  electrons,  thus 
the  term  self-sustaining.  The  voltage  at  which  a  spark  will  occur,  Vg,  is  a 
function  of  the  product  of  gas  pressure  and  gap  distance,  pD,  The  functional 
dependence  of  Vg  on  pD  is  described  by  a  Paschen  curve  which  is  peculiar 
to  the  discharged  gas.  The  Paschen  curve  for  hydrogen  is  shown  in  Fig.  2.  1 
and,  as  an  example,  it  can  be  seen  at  1  atm  and  a  gap  length  of  1  cm  an 
applied  field  over  10,000  volts  is  required  for  sparking. 

There  are  a  variety  of  types  of  self- sustained  (spark)  discharges  ; 
however,  only  two  characteristic  forms  will  be  considered  in  the  present 
discussion.  The  first  of  these  is  a  glow  discharge,  characterized  by  elec¬ 
tron  densities  of  ~  lO^^/^c  and  electron  temperatures  of  several  eV,  and 
the  second  is  a  thermalized  arc  discharge  characterized  by  electron 
densities  of  ~  lO^^/cc  and  an  electron  temperature  of  order  one  eV,  in 
equilibrium  with  the  gas  translational  temperature. 
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Fig.  2,1  Paschen  Curve  for  Hydrogen 
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These  three  types  of  discharges,  Townsend,  glow  and  thermalized 
arc,  have  been  examined  under  generic  conditions  for  discharge  initiation 
in  order  to  determine  which  type  is  best  suited  for  LSC  ignition.  This  is 
not  meant  to  imply  that  all  three  types  of  discharge  may  be  attainable  in  a 
specific  environment.  Indeed  Townsend  and  glow  discharges  are  more 
favorably  produced  at  sub-atmospheric  pressures;  however,  the  properties 
of  a  discharge  are  determined  not  only  by  pressure  but  also  by  numerous 
other  factors  such  as  voltage,  gap  size,  electrode  configuration,  and  the 
energy  content  of  the  relevant  electronics. 

The  scenario  considered  is  that  of  a  discharge  operating  perpendic¬ 
ular  to  the  gas  flow,  the  latter  being  co-axial  with  the  laser  beam.  The 
distance  between  anode  and  cathode  is  not  specified  but  must  be  larger  than 
the  laser  spot  size  since  laser  irradiation  of  the  electrode  elements  is  not 
desirable.  The  discharge  length  and  gas  velocity  have  arbitrarily  been 
chosen  to  be  1  cm  and  10^  cm/sec  respectively,  so  the  gas  flow  time  in 
the  discharge  region  is  1  msec.  Note  this  assumption  does  not  imply  that 
the  discharge  occurs  simultaneously  along  the  1  cm  length.  In  the  case  of 
a  spark  discharge  the  streamers  would  be  considerably  smaller  than  this. 
However,  once  a  spark  is  initiated,  it  will  attempt  to  travel  along  the  path 
of  least  resistance  and  thus  follow  the  flowing  gas  for  some  distance. 


Once  the  discharge  is  initiated  the  gas  in  the  discharge  region  will 
begin  to  absorb  the  laser  energy  via  neutral- electron  inverse  Bremsstrah- 
lung  (or  if  the  electron  density  is  high  enough  by  ion- electron  inverse 
Bremsstrahlung).  The  total  energy  absorbed  per  unit  volume  of  ionized 
gas  is  given  by 


AE  = 


/ 


k  I  dt 


(2.  1) 


where  I  is  laser  intensity  and  is  the  inverse  neutral  Bremsstrahlung 
absorption  coefficient  at  the  laser  frequency,  defined  by 

k^  =  aj_^  (T)  nn^  cm  ^  (2.2) 

where  n  is  total  gas  number  density,  ng  is  the  electron  number  density  and 
ttL  (^)  ^  function  of  temperature  (or,  more  appropriately,  electron 

temperature  in  non- equilibrium  situations)  and  frequency  which  is  specific 
to  the  absorbing  gas. 
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As  the  gas  absorbs  energy  it  will  heat,  causing  an  initial  pressure 
rise  and  subsequent  expansion  which  complicates  the  evaluation  of  relation¬ 
ship  (2.  1).  In  the  present  analysis  a  firm  upper  bound  on  the  energy 
absorbed  by  the  gas  has  been  determined  by  assuming  that  the  number 
density  of  the  gas  remains  constant  as  the  temperature  increases.  In  this 
case,  since  the  electron  density  of  the  discharge  is  specified,  k  will  be 
determined  by  the  form  of  a^,  (T).  In  the  cold  gas  the  dominant  neutral 
species  will  be  H2:  however,  as  the  gas  is  heated,  dissociation  will  occur 
and  hydrogen  atoms  will  become  the  dominant  neutral  partner  for  the 
absorption  process.  The  reduced  absorption  coefficients  kj^'  for  H  and  H2 
are  shown  in  Fig.  2.  2  as  a  function  of  temperature,  for  a  wavelength  of 
10.6  Hm.  These  predictions  have  been  developed  from  sources  described 
in  Reference  1.  It  is  to  be  noted  that  the  reduced  absorption  coefficients  for 
the  two  gases  are  quite  similar  and  that  they  scale  approximately  inversely 
with  the  temperature.  The  reduced  absorption  coefficient  may  be  related  to 
ttL  (T)  of  Eq.  (2.  2)  by  the  relationship 


(T)  kT  ,  (2.  3) 

where  k  is  the  Boltzmann  constant,  and  thus  for  these  gases  a (T)  is 
approximately  a  constant  for  \  =  10.6um.  Therefore,  from  Fig.  2.2,  k^ 
may  be  approximated  by 


k  «  3.  5  X 


10 


•37 


nn 


E’ 


cm 


(2.4) 
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Taking  n  =  2.  5  x  10  part/cc,  corresponding  to  p  =  1  atm  and  T  =  300°K, 
Eq.  (2.  1)  then  reduces  to 


AE 


10-2®  n^I, 


Joules/cc 


(2.5) 


Equation  (2.  5)  provides  an  approximate  estimate  of  the  amount  of 
energy  absorbed  in  one  msec  by  a  constant  density,  constant  electron  density 
plasma.  The  increase  in  temperature  of  the  gas  which  is  a  result  of  this 
absorption  may  be  deduced  by  consideration  of  the  thermodynamic  properties 
of  heated  hydrogen.  The  variation  with  temperature  of  the  enthalpy /unit 
mass  of  constant  pressure  hydrogen  plasmas  is  shown  in  Fig.  2.  3  for 
pressures  between  1  -  1000 atm.  These  predictions  were  taken  from  the 
analysis  of  Patch^.  The  enthalpy  variation  of  a  constant  density  hydrogen 
plasma  may  be  deduced  from  a  cross  interpolation  of  the  constant  pressure 
data.  Thus  for  a  given  amount  of  absorbed  energy,  as  predicted  by  Eq.  (2.  5), 
a  unique  plasma  temperature  may  be  determined. 
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Neutral  Inverse  Bremsstrahlung  Reduced 
Absorption  Coefficient  for  10,6  |j,m 


298  (  J  /  ^9 ) 


Squation  (2.  5)  has  been  evaluated  as  a  function  of  laser  intensity  for 
the  three  generic  discharges  discussed  earlier.  The  results  are  shown  in 
Fig.  2.  4.  The  final  hydrogen  plasma  temperature  which  would  be  attained 
for  a  given  amount  of  absorbed  energy  is  shown  on  the  right  hand  scale. 

Also  shown  as  tick  marks  on  the  curves  is  the  temperature  at  which  the  level 
of  thermal  ionization  exceeds  the  ionization  level  created  in  the  discharge. 
This  latter  temperature  is  a  useful  criterion  for  ignition  inasmuch  as  a  stable 
LSC  wave  would  exhibit  ionization  levels  higher  than  that  created  in  the 
discharge  (at  the  specified  total  densities).  It  can  be  seen  from  Fig.  2.4 
that  a  Townsend  discharge  will  not  produce  ignition  and  that  a  glow  discharge 
would  be  effective  at  best  only  for  intensities  >  10^  W/cm^.  Recalling  that 
the  results  of  Fig.  2.4  represent  a  strong  upper  bound  on  the  absorbed 
energy,  it  would  appear  that  the  thermalized  arc  discharge  is  best  suited 
for  the  successful  electric  discharge  initiation  of  LSC  waves  in  hydrogen. 

From  a  systems  viewpoint  it  is  desirous  to  define  the  minimum 
discharge  requirements  for  successful  ignition,  in  terms  of  electrical 
energy  content,  device  size  and  sophistication.  Unfortunately,  it  is  difficult 
to  provide  such  a  definition  by  theoretical  means.  The  production  of  stable 
thermalized  arcs  is  still  more  of  an  art  than  a  science.  One  standard  tech¬ 
nique  for  their  generation  is  through  use  of  an  over-voltaged  capacitor.  The 
evidence^  suggests  that  a  plasma  with  electron  densities  of  ~  10l5/cc  is 
forced  from  the  cathode  end  by  high  velocity  ionizing  waves  into  a  spark 
channel.  This  streamer  provides  a  finite  resistance  path  for  discharge  of 
the  capacitor.  As  the  capacitor  discharges,  the  gas  in  the  streamer  region 
will  be  heated  by  ohmic  power  dissipation.  Presently,  this  behavior  cannot 
be  predicted  a  priori  inasmuch  as  the  plasma  resistivity  is  a  function  of 
plasma  electron  density  and  streamer  size,  and  the  electron  density  in  turn 
is  connected  to  the  rate  of  ohmic  heating.  This  analysis  is  further  compli¬ 
cated  by  the  hydrodynamic  behavior  of  the  gas^.  The  initial  spark  can  lead 
to  formation  of  a  cylindrical  shock  wave  which  can  produce  further  ionization. 

In  any  event,  such  discharges  are  readily  produced  in  hydrogen  and 
thermal  equilibration  times  have  been  observed  to  be  tens  of  nanoseconds^* 
The  discharged  gas  is  low  density  (relative  to  pj^)  but  highly  ionized.  Indeed 
measurements  '  performed  in  hydrogen  in  near  atmospheric  discharges 
exhibit  electron  concentrations  near  the  anode  as  high  as  lO^^/cc  and  final 
temperatures  as  high  as  50  -  60,  000°K.  It  is  concluded  that  LSC  wave 
ignition  may  be  achieved  through  use  of  a  thermalized  arc  discharge;  how¬ 
ever  ,  it  is  felt  that  optimum  discharge  conditions  should  be  determined  by 
laboratory  experiments  rather  than  through  theoretical  analysis. 


-11- 


3.  PLASMA  INITIATION  BY  LASER-INDUCED  BREAKDOWN 


3.  1  INTRODUCTION 

Free  electrons  can  absorb  laser  radiation  during  collisions  with 
heavy  particles  such  as  ions,  atoms  or  molecules.  In  this  process,  known 
as  inverse  Bremsstrahlung,  the  role  of  the  heavy  particle  is  to  conserve 
momentum.  Little  of  the  energy  of  the  absorbed  photon  is  transferred  to 
the  heavy  particle.  Free  electrons  transfer  energy  among  themselves 
easily  during  collisions;  therefore  the  laser  energy  absorbed  by  the  elec¬ 
trons  is  rapidly  distributed  among  all  the  free  electrons.  The  electron 
distribution  can  then  be  characterized  by  an  electron  energy  e  roughly  equal 
to  two-thirds  the  average  electron  energy  just  as  kT  is  related  to  two-thirds 
the  average  energy  in  a  Maxwellian  distribution,  (The  factor  of  two-thirds 
relating  e  and  the  average  energy  is  exact  only  for  Maxwellian  distribution; 
nevertheless  the  characteristic  energy  e  is  a  useful  concept,  ) 

Energy  transfer  from  the  electrons  to  the  heavy  particles  is  less 
efficient  than  electron-electron  transfer;  hence  it  is  possible  for  the  gas 
temperature  to  lag  behind  the  electron  energy.  For  strong  fields,  electron 
heating  occurs  so  rapidly  that  the  few  free  electrons  reach  high  characteris¬ 
tic  energies,  that  is,  a  sizable  fraction  of  the  heavy  particle  ionization 
potential,  without  heating  the  gas  noticably.  At  such  energies  a  significant 
number  of  the  electron  -  heavy  particle  collisions  result  in  ionization,  A 
chain  reaction  then  occurs  with  both  the  original  electron  and  the  newly 
created  electron  absorbing  laser  energy  and  participating  in  further  ionizing 
interactions.  After  a  short  time  the  avalanche  created  by  electron  -  impact 
ionization  leads  to  a  sizable  concentration  of  high  energy  electrons  in  a 
otherwise  cold  gas,  i,  e,  ,  to  gas  breakdown. 

Experimentally,  as  laser  intensity  is  increased  from  an  initial  low 
level  an  abrupt  threshold  is  passed  above  which  this  breakdown  occurs.  ® 

This  threshold  intensity  can  be  estimated  by  balancing  two  important  physi¬ 
cal  mechanisms.  The  first  balance  determines  the  number  of  electrons 
which  remain  in  the  laser  -  heated  region.  The  rate  of  production  of  new 
electrons  must  eqtial  the  rate  at  which  they  are  lost.  The  second  balance 
determines  the  characteristic  energy  the  electrons  reach.  They  are  heated 
by  the  laser,  but  lose  heat  by  collisions  with  the  gas  molecules.  From  an 
analysis  of  there  two  balances,  an  estimate  of  the  intensity  for  laser  break¬ 
down  in  hydrogen  will  be  obtained. 
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3.  2  ELECTRON  PRODUCTION  AND  LOSS 


In  addition  to  production  of  electrons  through  electron  -  impact 
ionization,  electrons  may  be  freed  by  multiquanta  absorption,  the  process 
in  which  an  atom  or  molecule  absorbs  many  low  energy  photons  rather  than 
one  high  energy  photon.  Whenever  an  appreciable  number  of  electrons  is 
present  electron  -  impact  ionization  dominates.  Multiquanta  absorption  is 
important  only  as  a  potential  source  of  the  initial  electron  population 
required  for  electron  -  impact  ionization  to  dominate.  ^ 

In  molecular  gases  the  loss  mechanisms  generally  are  attachment  of 
an  electron  to  a  neutral  molecule,  recombination  of  an  electron  with  an 
ionized  species  (with  or  without  the  presence  of  a  third  particle)  and  electron 
diffusion  out  of  the  volume  influenced  by  the  laser.  Attachment  is  unimpor¬ 
tant  in  hydrogen.  The  recombination  rate  is  also  small,  particularly  during 
the  initial  stages  of  breakdown.  Diffusion  is  the  predominant  loss  process. 

The  net  rate  of  electron  production  in  hydrogen  can  thus  be  expressed 
as 


b  t 


1  E 


+  7- 


(3.  1) 


where  nj^  is  the  number  density  of  free  electrons,  is  the  ionization  rate 
per  electron  (the  concentration  of  neutral  molecules  is  absorbed  into  the 
definition  of  v.)  and  D  is  the  diffusion  coefficient. 

For  rough  estimates  the  diffusion  loss  mechanism  can  be  approxima¬ 
ted  by 


Dn 


E 


Ri 


(3.2) 


where  R^  is  the  radius  of  the  laser  beam.  As  long  as  the  electron  density 
is  low  the  electron  diffuses  freely.  At  high  concentrations  diffusion  proceeds 
at  the  slower  ambipolar  rate.  For  our  estimates  of  threshold  intensity  we 
will  use  the  free  diffusion  coeffient. 

Breakdown  requires  electron  production  to  exceed  diffusion  losses. 
The  threshold  electron  energy  for  breakdown  can  be  estimated  by  balancing 
losses  and  gains. 
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(3.3) 


Vi  (e)  =  (e) 


The  quantities  V^/nj^,  Dnj^  and  e,  where  nj^  is  the  number  density  of  free 
rnolecules,  are  known  as  functions  of  the  electric  field  per  particle,  E/n2^, 
from  experiments  with  static  electric  fields.  Despite  the  differ¬ 

ence  in  the  electron  distribution  function  for  DC  fields  as  compared  to  AC 
fields  calculations  have  indicated  that  the  ionization  rate  for  the  same  e  is 
very  similar.  We  have  thus  used  the  DC  experimental  radius  to  determine 
Vi/nM  and  Dnj^  as  functions  of  electron  energy  e.  In  Fig.  3.  1  we  have 
plotted  versus  e  the  diffusion  frequency  per  molecule  for  an  electron, 
Vd/^m*  for  various  values  of  pR^^,  (where  p  =  (3.  034E-17)  njy[  at  293K), 
and  the  ionization  frequency  per  molecule  for  an  electron,  v^/nj^.  The 
extremely  steep  increase  of  V^/nj^  with  e  indicates  why  the  threshold  is 
abrupt.  A  small  increase  in  e  above  the  temperature  needed  to  balance 
diffusion  losses  results  in  ionization  completely  dominating  diffusion.  For 
example  ionization  and  diffusion  balance  at  e  ~  2.  7  eV  when  p  =  1  torr- 
cm.  If  the  electron  temperature  reaches  2.  8  eV  for  the  same  pRj^,  we  find 
that 


— ^  -  — —  ~  2  X  10  cm^/sec  .  (3.4) 

M  M 

19  -3 

For  1  atmosphere  hydrogen  at  293K,  ~  2.  5  x  10  cm  so  that  solving 

the  electron  production  Eq.  (3.  1)  leads  to 


n  ~ 


n^  exp  (5.  0  X  10 


t) 


(3.5) 


Thus  if  the  breakdown  requirements  are  exceeded  by  even  a  small  amount 
the  gas  breaks  down  rapidly.  [Of  course  other  time  scales  are  also 
involved  -  the  time  to  heat  the  electrons  and  the  time  required  to  produce  a 
reasonable  initial  electron  population  if  one  is  not  present.] 

3.  3  ELECTRON  ENERGY  BALANCE 

The  minimum  laser  intensity  necessary  to  reach  the  breakeven 
electron  energy  is  determined  by  equating  the  electron  heating  by  the  laser 
and  the  power  lost  through  collisions  with  the  molecules.  The  electron 
heating  is  described  by 
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Fig.  3.1  Diffusion  Frequency  and  Ionization  Fr< 

in  Hydrogen  at  Various  Values  of  pR^ 
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dg 

dt 


(3.6) 


V.  e 

in 

where  e  and  m^,  are  the  charge  and  mass  of  the  electron,  respectively,  UJ 
is  the  angular  frequency  of  the  laser  radiation,  E  is  the  RMS  electric  field, 
momentum  transfer  collision  frequency  per  electron  and  ^he 
inelastic  collision  frequency  per  electron.  The  first  expression  on  the 
right  hand  side  represents  inverse  Bremsstrahlung  which  can  be  calculated 
from  collision  theory.  For  our  purposes,  however,  it  is  simpler  to  use 
the  expression  given  above  and  the  experimentally  determined  values  of 
Vm  and  Vin  versus  e  as  determined  from  DC  experiments. 

The  minimum  laser  intensity  I  necessary  to  reach  e  is  found  by 
setting  the  right  hand  side  of  Eq.  (3.6)  equal  to  zero. 

We  have 

I  =  cE^/4rr 


For  I  in  watts /cm 


2 


and  e  in  eV,  Eq.  (3.  7)  becomes 


I  =  1.  5  X  10 


2  2 
( UD  +  V  )  e  . 
m 


(3.8) 


It  is  now  possible  to  investigate  qualitatively  how  close  the  breakeven  inten¬ 
sity  is  to  the  true  threshold  intensity.  In  the  eV  range  increases 
extremely  rapidly  with  small  increases  in  G  where  as  is  almost  constant. 
Hence  the  second  term  of  Eq.  (3.  6)  is  negligible  throughout  most  of  the 
heating.  It  serves  only  to  stop  the  heating  abruptly  at  the  final  electron 
energy.  Because  of  the  sensitivity  of  both  the  production  rate  and  the  heat¬ 
ing  rate  to  small  changes  in  G,  the  two  balances,  Eqs.  (3.  3)  and  (3.6) 
should  lead  to  accurate  estimates  of  the  threshold  intensity.  There  is 
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greater  vincertainty  in  the  experimental  values  of  Vi^,  Dnj^  and 

e  then  in  the  identification  of  the  threshold  intensity  with  the  breakeven 
intensity. 

3.  4  BREAKDOWN  ESTIMATE 

For  a  given  beam  radius  and  pressure  p,  the  laser  intensity  to 
produce  breakdown  can  now  be  estimated.  Using  Fig.  3.  1,  the  energy  e 
at  the  crossing  point  of  the  and  Vj^/nj^  curves  is  found.  Using  this 

energy,  the  quantities  V^i  are  obtained  from  Refs.  10  or  11.  Then 

the  laser  intensity  I  is  calculated  for  the  desired  circular  frequency  (d 
from  (3.  8).  Using  this  procedure  the  threshold  intensity  necessary  to 
achieve  breakdown  in  pure  hydrogen  initially  at  293°  K  and  3  atmospheres 
pressure  is  plotted  versus  laser  beam  radius  for  a  10.6  (J.  laser  in  Fig.  3.2. 

The  threshold  intensity  is  seen  to  decrease  with  increasing  spot 
size  until  a  plateau  near  I  =  4.  5  x  10*^  W/cm^  is  approached.  The  decrease 
is  caused  by  the  lowering  of  the  diffusion  loss  rate  with  increasing  radius. 
The  plateau  represents  the  minimum  intensity  needed  to  raise  the  electron 
energy  high  enough  for  appreciable  ionization  to  occur.  The  intensity  of  a 
10  kW  laser  is  also  plotted  as  a  function  of  beam  radius.  Over  a  limited 
range  of  pressure  and  wavelength  the  breakdown  curve  scales  as  follows; 
for  a  given  laser  intensity  and  wavelength,  the  beam  radius  is  inversely 
proportional  to  pressure;  for  a  given  radius  and  pressure  the  threshold 
laser  intensity  is  proportional  to  the  wavelength  squared. 

Focussing  a  powerful  laser  unto  a  small  volume  will  produce  a  high 
concentration  of  free  electrons  within  a  small  volume.  The  effectiveness  of 
laser  induced  breakdown  in  initiating  an  LSC  wave  is  probably  not  limited  by 
the  electron  concentration  achievable  but  rather  by  the  volume  over  which 
such  a  concentration  can  be  sustained. 

3.  5  LASER-INDUCED  BREAKDOWN  WITH  PARTICULATES 

The  effect  of  particle  contaminants  on  the  breakdown  threshold  for 
air  is  a  subject  which  has  recently  received  considerable  attention.  It  has 
been  found  by  numerous  investigator s  14  -  16  that  the  intensity  threshold  in 
the  presence  of  particulates  may  be  lowered  as  much  as  three  orders  of 
magnitude  from  the  clean  air  threshold.  Although  there  is  no  data  for  the 
effects  of  particulates  upon  the  breakdown  threshold  in  hydrogen,  the  mech¬ 
anisms  by  which  particules  induce  breakdown  at  lower  intensities  in  air  will 
also  lead  to  lowering  of  the  threshold  in  hydrogen. 
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LASER  INTENSITY  I  (W/cm^) 


SPOT  RADIUS  R|j  (cm) 


Fig.  3.2  Laser  Intensity  Required  for  Breakdown  in  Hydrogen 
and  Laser  Intensity  Available  from  10  kW  Laser 
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A  summary  of  the  observed  breakdown  data  in  air  is  presented  in 
Fig,  3,  3.  The  air  breakdown  threshold  is  a  function  of  pulse  time,  beam 
diameter  and  particle  size.  The  lowest  observed  threshold  was  of  the 
order  of  a  megawatt  per  cm^;  however,  30iJm  particles  were  required  to 
initiate  the  breakdown.  In  addition,  it  was  observed^^  that  no  breakdown 
occurred  with  the  larger  particles  at  the  lowest  intensities  until  lOjisec 
into  the  pulse.  Shorter  pulses  do  not  yield  breakdown  with  contaminants 
xmless  the  flux  is  considerably  higher. 

The  physics  of  particle  induced  breakdown  is  not  well  understood. 

For  short  pulses  the  particles  evaporate  rapidly,  and  cascade  breakdown 
occurs  in  either  the  particle  vapor  or  the  air  surrounding  the  particle. 
Triplett  and  Boni^S  have  performed  a  theoretical  study  of  the  interaction 
of  laser  radiation  with  a  single  suspended  particle.  Their  analysis  applies 
only  when  the  laser  intensity  is  sufficiently  large  that  the  rapidly  evapora¬ 
ting  or  ’’exploding'*  particle  drives  a  strong  shock  wave  into  the  surrounding 
air.  Breakdown  then  occurs  via  a  non  equilibrium  cascade  process  (des¬ 
cribed  above)  in  either  the  vapor  or  shock  heated  air.  Calculations  indicate 
that  reductions  in  the  breakdown  threshold  to  slightly  greater  than  one  order 
of  magnitude  below  that  of  clean  air  are  feasible  by  this  mechanism.  It  is 
not  possible  to  predict  with  this  mechanism  the  reduction  in  the  threshold  by 
the  observed  several  orders  of  magnitude  at  longer  laser  pulses.  In  addi¬ 
tion,  an  ’’incubation”  time  or  delay  time  between  laser  turn-on  and  break¬ 
down  was  observed  with  the  long  laser  pulses^^  (  10  jJlsec).  Since  the  delay 

times  are  10  -  30  (isec,  the  breakdown  mechanism  is  nonequilibrium  cascade 
ionization  of  shock  heated  air  or  vapor.  It  appears  that  the  vapor  formed 
by  the  evaporating  particle  heats  in  local  thermodynamic  equilibrium  and 
transfers  its  energy  to  the  surrounding  air  by  conduction  and  radiative 
transport.  The  air  ultimately  heat  via  inverse  Bremsstrahlung  absorption 
until  the  formed  plasma  becomes  opaque  to  laser  radiation.  However,  at 
present  there  are  not  theoretical  models  capable  of  predicting  the  break¬ 
down  threshold  under  these  conditions. 

Therefore,  although  it  is  not  possible  to  predict  the  breakdown 
threshold  in  particle  seeded  hydrogen,  we  would  expect  this  to  be  a  viable 
technique  for  igniting  the  plasma  at  lower  intensities  than  predicted  by 
Fig.  3. 3. 
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Fig.  3.3  Measurements  of  Air  Breakdown  in  the 
Presence  of  Particulate  Matter 
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4.  0  PLASMA  INITIATION  BY  MOLECULAR  AND 
PARTICULATE  ABSORPTION 


The  two  methods  of  initiation  discussed  in  the  previous  two  sections, 
electric  discharge  and  laser-induced  breakdown,  provide  the  electrons  to 
initiate  inverse  Bremsstrahlung  absorption  by  a  nonequilibrium  process. 

A  third  alternative  is  to  provide  the  electrons  in  an  equilibrium  process  by 
heating  hydrogen  into  its  ionization  region.  For  pure  hydrogen,  that  means 
heating  to  T  ~  10,  000  K.  If  the  hydrogen  is  seeded  with  an  easily  ionizable 
material,  a  much  lower  temperature  would  suffice.  Cesium,  for  example, 
ionizes  around  3000  K. 

Such  heating  could  be  accomplished  if  one  could  find  molecules  or 
particles  which  would  absorb  lO.biim  laser  radiation,  become  hot  and 
transfer  their  heat  to  the  gas.  In  Ref.  1,  Sec.  3,  a  study  was  made  of 
possible  molecular  hbsorbers  for  10.  bjim  radiation.  It  was  concluded  that 
high  temperatures  could  not  be  reached  because  the  molecular  absorbers 
dissociate  (and  so  stop  absorbing)  at  a  few  thousand  degrees.  No  mole¬ 
cules  which  circumvent  this  limitation  were  identified  in  Ref.  1,  nor  have 
any  come  to  light  during  performance  of  the  work  reported  here.  However, 
if  seeded  hydrogen  is  considered,  then  the  gas  need  only  be  heated  to  a  few 
thousand  degrees  for  ionization  to  occur.  Molecular  absorption  might  be 
an  effective  means  of  accomplishing  this  heating. 

Particulate  absorption  was  also  studied  in  Ref.  1,  Sec.  4,  as  well 
as  in  Ref.  19.  Here  the  particles  absorb  laser  photons,  get  hot,  and  con¬ 
duct  heat  to  the  hydrogen.  Here  again  a  temperature  limitation  was  found. 
An  upper  bound  for  particle  vaporization  is  5000  K,  limiting  the  temperature 
which  can  be  reached  to  below  the  pure  hydrogen  ionization  region,  but 
within  the  ionization  region  of  seed  materials. 

There  was  noted  in  Refs.  1  and  19  some  possible  benefits  from  hot 
particles,  other  than  heating  the  gas.  The  vapor  produced  from  hot  parti¬ 
cles  may  itself  absorb  the  laser  radiation,  and  may  also  have  a  lower 
ionization  potential  than  the  hydrogen.  Furthermore,  small  particles  can 
emit  electrons  thermionic  ally  as  they  are  heated,  though  this  emission  is 
inhibited  by  electrostatic  effects  as  the  particle  charge  builds  up.  However, 
the  electrons  produced  may  be  sufficient  to  support  laser-induced  break¬ 
down. 


Nevertheless,  it  seems  at  the  present  time  that  molecular  or  parti¬ 
culate  absorption  is  not  a  favorable  scheme  for  initiating  an  unseeded 
hydrogen  plasma,  compared  to  electric  arc  or  laser-induced  breakdown. 
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5.  RADIATIVE  EMISSION  FROM  HYDROGEN 


5.  1  INTRODUCTION 

As  the  laser  heats  the  hydrogen  to  high  temperatures  the  radiation 
from  the  gas  increases  dramatically.  In  an  ionized  plasma,  radiation  con¬ 
tributes  importantly  as  an  energy  transfer  mechanism  within  the  plasma 
and  as  a  loss  mechanism  from  the  plasma.  The  radiative  transfer  within 
the  gas  affects  both  the  axial  and  radial  temperature  gradients.  Thus  a 
precise  description  of  the  physical  state  of  the  gas  in  the  rocket  core 
requires  that  the  radiation  and  fluid  mechanics  be  coupled.  It  is  beyond  the 
scope  of  this  report  to  deal  with  all  aspects  of  the  coupled  problem.  In 
keeping  with  approximations  to  be  made  in  the  fluid  mechanics,  such  as 
the  omission  of  the  radial  temperature  gradients,  we  neglect  radial  radia¬ 
tive  transfer  within  the  plasma.  Further-more,  the  axial  radiative  trans¬ 
fer  within  the  hot  gas  is  expected  to  be  small  compared  to  the  laser  heating; 
it  will  not  be  considered  here.  Axial  radiative  transport  from  the  hot 
plasma  to  the  cold  incoming  gas  may  be  important.  Its  effect  is  included  in 
the  LSC  wave  properties  discussed  in  the  next  section. 

The  only  plasma-produced  radiation  remaining  to  be  considered  is 
the  radial  flux  leaving  the  gas.  This  radiative  transport  enters  into 
rocket  engine  analysis  in  two  ways.  First,  the  power  lost  radially  must 
be  absorbed  somewhere  else  in  the  rocket,  either  in  the  walls  or  in  a  buffer 
gas.  Second,  the  radial  radiation  which  escapes  the  plasma  must  be 
accounted  for  in  the  fluid  mechanical  equations  as  a  loss.  At  sufficiently 
high  temperatures,  in  fact,  it  is  the  major  loss  mechanism  and  effectively 
limits  the  maximum  temperature  attainable. 

The  approach  we  follow  here  is  first  to  study  the  local  radiation 
properties  of  the  plasma,  and  then  to  develop,  using  the  local  properties 
as  guides,  a  reasonable  approximation  to  the  coupled  radiating  gas  dynamic 
system. 

5.  2  LOCAL  SPECTRAL  RADIATION  PROPERTIES 

The  spectral  radiation  properties  of  a  hydrogen  plasma  are  deter¬ 
mined  by  the  local  state  of  the  system.  We  assume  that  the  plasma  is  in 
local  thermodynamic  equilibrium,  so  that  temperature  and  pressure  define 
the  state  of  the  gas.  The  concentrations  of  the  various  species  is  determined 
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through  the  equilibrium  conditions.  The  temperatures  and  pressures  of 
interest  are  those  expected  in  a  rocket  heated  by  laser  absorption;  that  is, 
temperatures  between  10,  000  K  and  20,  000  K  and  pressures  between  one 
and  ten  atmospheres.  In  the  interest  of  simplicity  we  need  only  include 
effects  which  are  expected  to  contribute  10%  or  more  of  the  total  radiation. 

The  basic  ingredients  for  determining  local  radiative  properties  are 
the  cross  sections  for  the  possible  interactions  and  the  number  density  of 
the  various  species.  The  use  of  local  thermodynamic  equilibrium  assumes 
that  collisions  among  the  various  species  of  the  gas  are  frequent  enough 
to  ensure  an  equilibrium  distribution  of  velocity  and  of  species,  corres¬ 
ponding  to  the  local  temperature  and  pressure.  With  this  assumption  and 
the  principle  of  detailed  balance  the  radiative  transfer  equation  can  be 
written  in  a  simplified  form  involving  only  one  function  depending  upon  the 
concentration  of  species  and  the  interaction  cross-section.  The  radiative 
transfer  equation  has  the  following  form:^® 


(x,  'n^)  =  k'  (T,  p)  [B  (T)  -  I  (x,  ^)  ]  (5.1) 

U)  UJ  UJ  ou 

where  au  =  hp  v  is  the  photon  energy;  is  the  spectral  radiation  intensity 
which  depends  on  UJ,  x,  and  the  direction  specified  by  the  unit  vector  n; 
and  (T)  is  the  equilibriunn  spectral  radiation  intensity  given  by 


‘  B  (T)  = 
0) 


ii) 


2  T-  3  /  U)/kT  3 
c  TT  (e  -1)4tt 


(5.2) 


where  c  is  the  velocity  of  light,  lx  is  the  Planck  constant  divided  by  2TT  and 
k  is  the  Boltzmann  constant.  The  effective  absorption  coefficient  is  related 
to  k^  (T,p),  the  absoption  coefficient  uncorrected  for  reabsorption,  by 


-U)/kT 

k^  (T,  p)  =  (1  -  e  )  k^  (T,  p)  .  (5.  3) 

The  absorption  coefficient  k^  (T,p)  is,  in  turn,  easily  expressed  as  a  sum 
of  products  of  absorption  cross  sections  and  species  concentrations.  The 
temperature  and  pressure  dependence  of  k^j  enter  mainly  through  the  con¬ 
centration  of  species.  However,  secondary  temperature  variations  occur 
whenever  the  cross  section  must  be  averaged  over  the  relative  flxix  of  some 
of  the  involved  species,  as  occurs  in  inverse  Bremsstrahlung. 
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The  first  step  in  determining  (T,p)  is  to  detail  which  interactions 
are  expected  to  be  important  for  the  range  of  temperatures  and  pressures 
expected.  The  interactions  to  be  considered,  in  turn,  depend  upon  the 
species  which  are  abundantly  present.  When  radiation  is  an  important  loss 
mechanism  the  concentration  of  molecules  is  several  orders  of  magnitude 
lower  than  the  atomic  concentrations.  We  therefore  consider  a  plasma  in 
which  only  atoms,  excited  atoms,  ions  (protons)  and  electrons  are  present. 
Their  relative  abundance  is  determined  through  the  Saha  equation 


v3/2 


n. 


^el  I  ^el  E  I  ^ 

‘elA  \2TTTi‘ 


3/2 


(5.4) 


where  nj  and  are  respectively  the  concentration  of  electrons,  ions 
and  hydrogen  atoms,  ^ell  ^elA  electronic  partition 

functions  of  the  electron,  ion  and  the  neutral  atom  and  0j  is  the  ionization 
portntial  expressed  in  temperature  units,  (In  more  accurate  calculations 
than  are  necessary  here,  0j  would  be  replaced  by  0j  -  A  0j  where  A  0j  is 
the  effective  lowering  of  the  ionization  potential  corresponding  to  termina¬ 
ting  the  sum  in  the  electronic  partition  function  of  the  neutral  atom  after  a 
finite  number  of  excited  states). 

The  transitions  which  lead  to  emission  can  now  be  listed.  They  are: 

1)  capture  of  a  free  electron  by  an  ion 

2)  Bremsstrahlung  involving  electron-ion  collisions 

3)  Bremsstrahlung  involving  electron-atom  collisions 

4)  radiative  decay  of  an  excited  atom. 

5.3  CONTINUUM  TRANSITIONS 

The  spectral  properties  of  the  first  three  types  of  transitions,  the 
transitions  involving  the  electron  continuum,  can  be  estimated  by  semi- 
classical  formulas.  The  quanturn  corrections  to  these  formulas  need  only 
be  included  when  the  quantum  correction  is  large  and  the  radiation  emitted 
is  a  sizeable  portion  of  the  total  radiation  produced. 

From  Kramer's  semi-classical  analysis  of  electron-proton  Brems¬ 
strahlung  the  absorption  coefficient  is  given  by^^ 
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where  e  is  the  charge  of  the  electron  and  the  pressure  dependence  enters 
only  through  the  number  densities. 


Similar  semi- classical  formulas  can  be  derived  for  bound-free 

absorption  wherein  an  electron  of  energy  in  the  nth  quantum  state  of 

an  atom  absorbs  a  photon  of  energy  UJ  resulting  in  a  free  electron  kinetic 

energy  ou  +  E  and  a  free  ion.  The  semi-classical  cross  section  for  absorp- 
,  ,  oo  n 

tion  18^^ 
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The  total  bound  free  absorption  coefficient  for  radiation  of  given  cu  is  then 


n  S:  n 
(1) 


a 

tt)  >  n 


(5.  7) 


where  is  the  number  density  of  atoms  in  the  nth  quantum  state,  and  where 
nyj*  is  the  lowest  quantum  state  which  can  be  ionized  by  a  photon  of  energy  au; 
that  is 


k  9j, 

(n  ) 
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n  -  1) 
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(5.8) 


The  Bremsstrahlung  due  to  ion-electron  collisions  Eq.  (5.  )  can  be 
re-expressed  in  a  form  similar  to  Eq.  (5.  7)  by  using  the  Saha  Eq.  (5.  4) 
with  Qeil  =  1  to  convert  the  product  n^  nj  into  the  concentration  of  ground 
state  atoms,  n^Q  and  by  similarly  expressing  n^^  in  terms  of  n^Q: 
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(5.10) 
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The  sum  of  the  absorption  coefficients  for  Bremsstrahlung  and  free- 
bound  transitions  is 


k  (T,  p)  = 
uu 


o  10 
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where  we  define  - ,  and  we  have  used  in  (5.  5)  the  expression  for  the 

n  T 

ionization  potential  of  hydrogen  in  terms  of  fundamental  constants. 


4  ,  r-  2 

=  e  k  h 


(5.12) 


to  get  the  first  term  in  brackets  in  (5.  11). 

.  -3  - 1 

If  uj  is  given  in  eV  and  in  cm  ,  then  k^^  in  cm  is 
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(If  the  lowering  of  the  ionization  potential  were  important,  Eq^  (5.  11)  would 

,  AGt/T  , 

have  1/2  Xj  replaced  by  e  /2X2). 

The  emission  coefficient  (T,p)  is  defined  as 


-Uj/kT 

=  k  (T,  p)  (1  -  e  )  B 

ou  ou 


(5.13) 
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The  semi-classical  approximation  to  (T,p)  found  from  Eqs.  (5.  11) 
and  (5.  13)  is  shown  in  Fig.  5.  1  for  p  =  3  atm  and  T  =  14,  000  K.  The  total 
continuum  power  emitted  per  unit  volume,  neglecting  reabsorption  effects, 
is  found  by  integrating  ©m  over  solid  angle  and  over  the  photon  spectrum. 
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(5.14) 


In  Fig.  5.  2,  (T,  p)  is  plotted  versus  temperature  for  a  hydrogen 

plasma  at  3  atm. 

If,  on  the  other  hand,  the  vacuum  ultraviolet  part  of  the  spectrum 
(tu  ^  13.  6  eV)  is  effectively  reabsorbed,  the  continuum  power  lost  in  the 
semi -classical  approximation  is 


13.6 


which  is  also  shown  in  Fig.  5.2. 

The  relative  magnitude  of  boimd-free  emission  and  Bremsstrahlung 
can  now  be  determined.  The  contribution  of  ion- electron  Bremsstrahlung  to 
the  emitted  power,  as  compared  to  the  bound-free  emission,  is  estimated  by 
comparing  1/xj  to  2.404  when  no  reabsorption  occurs  and  to  0.404  when  the 
vacuum  ultraviolet  is  strongly  reabsorbed.  For  the  temperatures  of  interest 
we  have  0.  07  <  <  0.  14.  Thus  Bremsstrahlung  is  not  expected  to  be  an 

important  factor  unless  the  bound  free  emission  is  inhibited  by  reabsorption. 

The  other  continuum  contribution,  Bremsstrahlung  associated  with 
electron-neutral  collisions,  is  much  smaller  than  electron-ion  Bremsstrah¬ 
lung  except  when  atoms  are  at  least  one  hundred  times  more  abundant  than 
ions.  Only  at  low  temperatures,  such  as  10,  000 K,  and  high  pressures 
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does  the  neutral  concentration  overwhelm  the  ion  concentration  sufficiently 
to  make  neutral  Bremsstrahlung  as  important  as  electron-ion  Bremsstrah- 
lung,  which  itself  is  only  a  small  fraction  of  the  total  radiation.  Moreover, 
as  can  be  seen  from  Fig.  5.  2,  plasma  radiation  is  very  small  under  these 
conditions.  For  our  purposes,  then,  the  contribution  of  neutral  Bremsstrah¬ 
lung  can  be  safely  omitted  as  an  important  radiation  loss  mechanism. 

The  expressions  quoted  above  are  semi-classical  in  origin.  To 
include  the  quantum  effects  the  cross-sections  must  be  multiplied  by  a 
correction  factor  known  as  the  Gaunt  factor.  For  free-free  transitions  at 
the  temperatures  and  pressures  imder  consideration  the  Gaunt  factor 
averaged  over  the  range  of  ou  is  approximately  1.4.  (The  precise  value 
varies  with  temperature  and  pressure. )  Since  this  is  a  correction  to  an 
already  small  contribution  it  is  not  included  in  further  computations. 

The  bound-free  transitions  also  have  Gaunt  factors  associated  with 
them.  The  most  important  quantum  correction,  about  20%,  concerns  the 
transition  to  the  ground  state.  Rather  than  using  a  Gaunt  factor,  however, 
it  is  more  convenient  in  this  case  to  begin  with  the  correct  quantum  cross 
section  for  the  region  of  u)  where  the  bulk  of  the  emission  occurs.  From 
Bethe  and  Salpeter24  ^e  have  to  replace  the  cross-section  of  (5.  6)  when 
n  =  1  by 
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(5.16) 


We  then  separate  off  the  continuum  n  =  1  term  of  the  sum  in  (5.  11) 
and  find  the  absorption  coefficient  to  be 
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(5.17) 


The  last  term  is,  of  course,  only  to  be  used  when  n^j  <  1.  This  quantum 
correction  is  the  only  one  we  include  in  the  calculations  in  Subsection  5.6. 
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5.4  ATOMIC  TRANSITIONS 


The  final  radiation  process  to  be  considered  is  the  transition  between 
two  atomic  bound  states.  The  total  power  emitted  in  a  spectral  line  corres¬ 
ponding  to  a  transition  from  an  upper  state  of  quantum  number  n'  to  a  lower 
state  designated  by  quantum  number  n  is  easily  computed  from  the  well 
known  f„„'  number  for  atomic  transitions. 


?  c 

The  f  n\imbers  of  the  most  important  lines^^  are  given  in  Table  5.1. 
The  f  number  of  the  higher  lines  (n'  >  8)  in  the  first  three  series  are  given 
by 
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(5.18) 


where  Aj  =  1.6  for  the  Lyman  series,  (n  =  1),  A2  ~  3.  7  for  the  Balmer 
series  (n  =  2),  and  A^  =  5.3  for  the  Paschen  series  (n  =  3).  The  f^n* 
numbers  for  transitions  involving  n>  3  can  be  approximated  as 


nn' 


(5.19) 


From  the  number  the  cross  section  integrated  over  photon  energy  is 
found  to  be  ^  ^ 

a  t  -  -  *  * 

^nn  c  im  , 

and  the  power  emitted  per  unit  solid  angle  is  given  by 


/ 


e  du)  =  (1  -  e 

ou 


2  -“in''''’' 

>  ,™°nn"’AG<“  ^  *' 

nn' 


line  oj  » 
nn 


where  uojj  =  k  0j 
be  rewritten  as: 


By  substituting  for  ,  (T),  Eq.  (5.20)  can 

nn 
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TABLE  5.  1 


f  NUMBERS  FOR  HYDROGEN  FOR  TRANSITIONS 
FROM  LOWER  STATE  n  TO  UPPER  STATE  n' 
FROM  REFERENCE  25 


n _ _ ^nn' 

1  2  .4162 

"  3  .0791 

”  4  .0290 

"  5  .0139 

"  6  .0078 

"  7  .0048 

"  8  .0032 

2  3  .641 

"  4  .1 195 

"  5  .0444 

"  6  .02235 

"  7  .0124 

"  8  .0082 

3  4  .8413 

”  5  .150 

”  6  .0554 

"  7  .0271 

"  8  .0162 
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The  power  emitted  (4tt  times  Eq.  (5.  21))  in  several  of  the  low-lying 
lines  dwarfs  the  total  continuum  contribution.  In  Table  5.  2  we  illustrate 
this  feature  for  hydrogen  at  14,  000  K  and  3  atm.  However  reabsorption, 
which  will  be  discussed  in  the  next  subsection,  reduces  the  magnitude  of 
the  losses  due  to  line  radiation.  In  order  to  adeq\iately  analyze  line  trans¬ 
port  it  is  essential  to  know  not  only  the  total  energy  radiated  but  also  the 
detailed  profiles  of  the  strong  lines.  The  profile  is  modified  by  the 
influence  of  other  particles.  In  the  presence  of  atoms,  ions  and  electrons 
the  shape  of  the  line  is  broadened  from  the  narrow  Lorentz  characteristic 
of  the  natural  line  shape;  however  the  change  in  total  power  emitted  is 
negligible.  In  an  ionized  plasma  such  as  in  the  hot  core  of  the  rocket  the 
broadening  is  caused  mainly  by  the  charged  particles  and  is  commonly 
termed  Stark  broadening.  Hydrogen  lines  are  affected  more  strongly  by 
Stark  broadening  than  lines  of  other  atoms,  so  that  in  ionized  hydrogen 
other  broadening  mechanisms  such  as  Doppler  broadening  are  unimportant. 

The  profiles  of  hydrogen  lines  subject  to  Stark  broadening  serve  as 
important  diagnostic  tools;  as  a  result  they  have  been  thoroughly  studied. 
Elaborate  calculations  of  the  profiles  have  been  performed  and  the  profiles 
carefully  tabulated.  For  our  purposes,  however,  only  the  far  wings  of 
the  Lyman  lines  are  important.  In  the  asymptotic  far  wings  the  profile  S 
is  adequately  described  by  2  times  the  Holtsmark  profile; 

S  (a)  =  2  C  ,/  la  I  (5.  22) 

where  C^n*  ^  constant  and  a  is  the  reduced  wavelength  defined  in  terms 
of  the  line  center  energy  tu^,  the  line  center  wavelength  and  the  normal 
field  strength 

2/3 

F  =e(4TTn  /3)  ,  a  =  (u)  -  («„  )  X  /u)  F 

o  i  ^  o  ^ 
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TABLE  5.2 


POWER  EMITTED  BY  LINES  IN  TRANSPARENT  APPROXIMATION. 
HYDROGEN  AT  3  ATMOSPHERES  AND  T  =  14000  K 


Power  Emitted 


Line 

n' 

n 

fwatts/cm^l 

La 

2 

1 

7.0  X  10® 

Lp 

3 

1 

4.62  X  10^ 

Ly 

4 

1 

2.4  X  lo"^ 

L  6 

5 

1 

4.59  X  10® 

6 

1 

2337 

7 

1 

1350 

8 

1 

866 

9  ^  00 

1 

2463 

Ha 

3 

2 

5705 

m 

4 

2 

1507 

Hy 

5 

2 

616 

6 

2 

314 

4 

3 

418 

5 

3 

176 

6 

3 

94.3 

5  00 

4 

80.4 

Continuum 

2. 6  X  10^ 
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The  profile  is  normalized  by  requiring 


00 

J's  (a)  da  =  1  (5.23) 

-00 

Thus  in  our  spectral  plots  of  the  emission  coefficient  vs  uj  we  have 


e  =  ^1 
OJ  du,  1 

S  (a)  /  e^du) 

(5.24) 

line 

^0 

■  S  (a)  /*  e  du) 

J  OJ 

(5.25) 

U^o  ^0 

line 

The  factor  of  two  in  Eq,  (5,  22)  is  required  in  order  to  include  the  electron 
broadening  effects  as  well  as  the  ion  broadening.  The  values  of  for 

the  first  three  Lyman  lines^7  are  given  in  Table  5,  3,  Typical  line 
broadened  profiles  are  illustrated  in  Fig,  5,  3,  Only  the  first  four  lines 
are  shown  isolated.  The  higher-lying  lines  are  indistinguishable.  They 
merge  to  form  an  extension  of  the  bound-free  continuum.  When  necessary 
more  detailed  line  profiles  can  be  found  in  Griem,  26 

5,5  SIMPLE  MODEL  OF  RADIATIVE  LOSSES 

Using  the  spectral  effective  absorption  coefficients  found  in  the  last 
subsection  we  can  now  proceed  to  model  the  radiation  from  the  hot  plasma 
core. 


In  principle  we  wish  to  find  the  spectral  intensity  I^  (x,  n)  at  each 
point  X  and  direction  "n"  in  the  core  by  solving  the  equation. 


'n-Vl  (x,  =  k'  (T,  p)  [B  (T)-I  (x, (5.26) 

U)  (1)  OU  UD 

subject  to  the  boundary  condition  that  at  the  outer  surface  of  the  plasma  the 
radiant  intensity  directed  into  the  plasma  be  zero.  Rather  than  attempt 
such  an  ambitious  task  here,  we  are  satisfied  to  develop  a  model  in  which 
the  plasma  radiation  can  be  treated  as  a  local  loss  term  in  the  gas  dynamic 
equations  and  in  which  the  total  flux  of  radiation  at  the  edges  can  be  estimated. 
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TABLE  5.3 


COEFFICIENTS  FOR  ASYMPTOTIC 


HYDROGEN  HOLTSMARK  PROFILES 
FROM  REFERENCE  27 


nn 


I 


5/2 


n 


n 


C  , 
nn 

angstrom 


cgs  field  strength  unit 


1 


3/2 


1  2  3.37  X  10 

1  3  1.80  X  10 

1  4  3.  99  X  10 
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cm^  eVsterad 


aj(eV) 

Fig.  5.3  Stark  Broadened  Profiles  of  Lyman  Series  in  Hydrogen 

and  the  Black  Body  Limits  for  the  Two  Choices  of  Radius 
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To  this  end  we  consider  tiie  radiation  lost  per  unit  volume  at  temperature 
T  and  pressure  p  from  the  plasma  to  be  equivalent  to  the  average  radia¬ 
tion  loss  per  unit  volume  from  an  infinitely  long  cylinder  of  uniform 
temperature  T,  pressure  p  and  radius  r.  In  this  way  the  major  qualita¬ 
tive  effects  of  reabsorption  are  estimated  without  including  the  effects  of 
temperature  variations. 


The  power  emitted  per  unit  length  at  a  given  u)  is 
P  =  2t?  B  (T)  R  F  [  k'  R] 

U)  (1) 

where  the  function  F  is  defined  as 


(5.27) 


Tt/2 

F^k^'R)=2y*  de  cos^e 


(5.28) 


The  function  Ij  and  Lj  are  respectively  the  modified  Bessel  function  of 
order  one  and  the  modified  Struve  function  of  order  one.  The  function 
F  (k(jy  R)  gives  the  ratio  of  actual  power  emitted  to  the  power  emitted  by 
a  blackbody  of  the  same  radius.  It  is  plotted  in  Fig.  5.4. 

Two  extreme  limits  are  important.  When  k^^  is  sufficiently  small 
so  that  k^  R  «  1  is  satisfied,  reabsorption  effects  are  unimportant. 

The  total  emission  per  unit  length  is  then  equivalent  to  the  emission  per 
unit  volume  multiplied  by  the  area; 


'’mlvolame)  =  ^  \  (5.29) 

t  * 

The  other  extreme  is  large  k^  .  When  k^  R  »  1  reabsorption  dominates 
and  the  cylinder  looks  like  a  blackbody  surface  radiator: 

,  ^0,  (surface)  =  =  (5-30) 

For  comparison  these  two  extreme  limits  are  also  plotted  in  Fig.  5.4  with 
the  factor  of  2tt^RB^  divided  out.  Calculations  using  Eq,  (5,25)  to 
weight  the  various  parts  of  the  spectrum  can  be  performed  numerically; 
however,  first  estimates  of  the  effect  of  radiation  can  be  made  by  resorting 
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Black  Approximation 


Ratio  of  Actual  Power  Lost  to  Black  Body  Limit  on  Power 
Lost  for  an  Infinite  Cylinder.  The  Transparent  and  Black 
Body  Approximations  Are  Also  Shown. 


to  an  even  stronger  approximation.  The  power  lost  per  unit  length  will 

be  approximated  as  the  smaller  of  P  ,  ,  ,  and  P 

(w  (volume)  tu  (surface) 

Using  Eqs.  (5.29)  and  (5.  30)  we  find  that  the  plasma  is  a  trans¬ 
parent  volume  radiator  for  all  u;  such  that  2  R  ^  1  is  true,  or 
equivalently,  such  that 


e 

U) 


(T,  p)  ^ 


2  R 


(5.31) 


is  satisfied.  The  transparent  radiation  losses  can  be  immediately  in¬ 
corporated  into  the  fluid  dynamic  equation  as  a  loss  term.  The  divergence 
of  the  plasma  flux  at  any  energy  ou  is  simply 


S  =  4  n  e  (T,  p)  .  (5. 32) 

If  the  opposite  inequality  is  valid,  i.e. 

B  (T) 

2  R  Si  (T,  P)  ,  (5. 33) 

the  plasma  acts  as  a  blackbody  radiator  with  effective  volume  radiation 

2rrB  (T)/R 

U) 

It  is  not  clear  whether  or  not  this  loss  should  be  included  in  the  fluid 
mechanics  as  an  effective  volume  loss.  Even  though  it  is  a  loss  from 
the  material,  the  loss  occurs  mainly  at  the  surface;  particularly  when 
2  R  »  1  is  satisfied.  Then  most  of  the  plasma  in  our  model  is  un¬ 
affected  since  the  temperature  is  assumed  to  be  held  fixed  throughout 
the  radius.  In  a  more  realistic  approximation,  radial  temperature 
gradients  must  exist  which  then  lead  to  heat  transfer  from  the  hotter 
interior  to  the  cooler  surface.  The  heat  transfer  is  accomplished  through 
both  thermal  conduction  and  radiative  transfer.  In  this  manner  the  surface 
losses  do  influence  the  interior.  However  a  realistic  radial  profile  may 
also  result  in  photon  trapping  especially  for  regions  of  the  spectrum 
where  k^^j  is  very  large.  Without  additional  information  concerning  the 
radial  temperature  profile  it  is  difficult  to  assess  accurately  the  role 
of  the  surface  radiation  in  fluid  equations,  especially  when  the  equations 
do  not  admit  radial  temperature  gradients.  This  ambiguity  suggests  that 
two  predictions  of  radiative  losses  should  be  made.  For  an  upper  limit  on 
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the  flux  which  leaves  the  plasma,  both  the  blackbody  and  transparent 
radiation  are  included.  For  a  reasonable  lower  limit  on  radiative  loss 
only  the  transparent  radiation  and  any  borderline  blackbody  radiation 
(say  2  k^'  R  <  3)  are  to  be  included. 


5,6  APPLYING  THE  MODEL 

Our  criterion  is  illustrated  in  Figs,  5,  1  and  5,  3  where,  in  addi¬ 
tion  to  the  contributions  to  the  emission  coefficient,  we  have  plotted  the 
function  (T)/{2R)  for  two  values  of  R  of  interest.  In  Fig,  5,  1  where 
the  classical  approximation  to  the  continuum  radiation  is  shown,  the  low 
energy  region  of  the  spectrum  is  seen  to  lie  below  the  curves  Bjjj(T)/2R 
for  both  radii.  Therefore  the  transparent  requirement  is  satisfied. 

The  vacuum  ultraviolet  portion  (uj  s  13,  6  eV)  lies  slightly  below 
Byj  (T)/2R  for  R  =  0,  05  mm  but  well  above  for  R  =  1  cm.  Hence  the 
VUV  is  barely  transparent  for  the  small  radius  and  strongly  black  for  the 
large  radius.  Four  distinct  Lyman  lines  are  shown  in  Fig.  5,  3.  The 
three  lowest  lines  are  strongly  black  at  their  peak  but  they  possess  wide 
transparent  wings.  The  higher  lines  tend  to  mimic  an  extension  of  the 
continuum.  In  our  calculations  of  power  lost  we  have  lumped  the  higher 
Lyman  lines,  L5  and  above,  with  the  true  continuum  insofar  as  deter¬ 
mining  whether  they  are  black,  marginally  black,  or  transparent.  To 
determine  the  power  lost  in  the  transparent  limit  for  the  region  of 
(U  ^  13,6  eV  we  have  used  the  quantum-mechanical  expression  (5,  17) 
for  the  continuum,  and  the  sum  of  line  strengths  for  n's  4  for  the  high 
lying  lines.  For  the  three  isolated  lines  the  centers  have  been  treated 
as  black.  Fortunately  the  transition  between  the  black  center  and  the 
transparent  wing  occurs  far  enough  from  the  line  center  that  the 
asymptotic  Holtsmark  profile  is  adequate  to  describe  the  line.  The 
transition  point  Aou  where  black  turns  to  transparent  can  be  determined 
by  finding  AU)  such  that  em  as  found  from  Eq.  (5.25)  satisfies  the  equality 
in  Eq.  (5.31). 


The  total  power  per  steradian  radiated  by  the  wings  is  thus  found 
by  integration  to  be 

c  .  r 


P  -  — 

line  wings  3 


-T',  I  e  du) 

3/2  J  I  00 

(arr.)  linenn 


(5.34) 


where  a-j-  is  the  reduced  wavelength  corresponding  to  the  transition  from 
black  to  transparent. 


The  only  other  line  which  is  strongly  reabsorbed  in  our  calcula¬ 
tions  is  the  first  Balmer  line  H^.  The  self  absorption  is  important 


-44- 


only  at  high  pressures  and  large  radii.  The  asymptotic  formulas  are 
not  sufficient  to  estimate  the  reabsorption,  thus  the  detailed  profiles 
of  Griem26  must  be  used. 

The  total  power  emitted  per  unit  volume,  including  all  transparent 
radiation  and  an  effective  vacuum  ultraviolet  radiation  whenever  it  is 
marginally  black,  is  illustrated  in  Figs.  5.  5  and  5.  6  as  a  function  of 
temperature  for  two  values  of  pressure  and  radius.*  In  the  numerical 
fluid  mechs-nical  calculations  of  Section  7,  the  radiation  loss  terms  are 
adequately  approximated  as  constant  factors  multiplying  one  of  the  semi- 
classical  formulas  (5.  14)  or  (5.  15).  These  approximations  are  also  in¬ 
dicated  in  Figs.  5.  5  and  5.6.  All  the  formulas  are  of  the  following 
form; 

P^=8.6x  l0-’  C^)  ,5.35) 

i  I  cm 

-3 

where  is  in  cm  .  The  constants  Cj^  and  C2  used  for  various 

pressures  and  radii  are  converted  to  a  different  form,  in  mks  units,  in 
Eq.  (7.38),  and  then  listed  in  Table  7,3. 

The  calculations  listed  above  are  reasonable  estimates  of  the 
effective  loss  term  per  unit  volume  from  the  bulk  of  the  hot  plasma. 
Radiation  which  is  subject  to  strong  reabsorption  has  not  been  included 
in  these  estimates.  The  blackbody  surface  contribution  can  be  easily 
included  by  integrating  u  (T)  over  the  region  of  the  spectrum  which 
satisfies  our  criterion  for  blackness.  In  order  to  compare  the  magnitude 
of  the  black  and  transparent  losses,  either  the  blackbody  losses  must  be 
rewritten  as  effective  volume  losses,  or  the  transparent  losses  must  be 
transformed  into  an  effective  surface  flux.  In  Table  5,4  we  have  com¬ 
pared  the  black  and  transparent  effective  surface  fluxes  for  two  situations 
of  interest. 

The  transparent  radiation  accounts  for  the  major  portion  of  the 
total  power  lost.  However  the  blackbody  contribution  is  still  a  significant 
fraction  of  the  total  loss.  Furthermore  the  fraction  of  the  total  emission 
which  is  blackbody  radiation  generally  increases  as  the  radius  becomes 


Radiation  which  is  black  is  not  included  because  of  uncertainty  in  the 
radial  temperature  gradients,  as  discussed  below  Eq,  (5,33). 
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T(IO^K) 

Power  Radiated  from  Hydrogen  (3  atm,  R  =  1  Showing 
Numerical  Results  and  the  Analytic  Approximation 


Fig.  5.5 


(  W  /cm^) 


Fig.  5,6 


Power  Radiated  from  Hydrogen  (10  atm,  R  =  .05  cm)  Showing 
Numerical  Results  and  the  Analytic  Approximation 


TABLE  5.4 

EFFECTIVE  SURFACE  RADIATION  LOSSES 


LASER 

10  kW 

5  MW 

Pressure 

3  atm 

3  atm 

Temperature 

14000  K 

14000  K 

Radius 

,056  cm 

1.25  cm 

Transparent 

2 

1579  watts /cm 

2 

6507  watts /cm 

Black 

2 

462  watts /cm 

2 

3356  watts /cm 

Total 

2 

2041  watts /cm 

2 

9863  watts /cm 
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larger.  For  accurate  results  at  larger  radii  it  is  necessary  to  include 
the  blackbody  radiation  properly,  but  this  is  difficult  unless  one  has 
knowledge  of  the  radial  temperature  profile,  as  discussed  below  Eq.  (5.  33). 

5,7  SUMMARY 

Having  analyzed  the  gross  behavior  of  the  spectral  emission  co¬ 
efficient  as  a  function  of  ou,  T  and  p,  we  have  obtained  estimates  of  the 
volume  power  loss  due  to  transparent  radiation,  and  of  the  black  body 
surface  radiation.  These  estimates  are  based  upon  a  simple  model  - 
one  in  which  the  radiation  lost  per  unit  volume  from  the  core  at  given  T 
and  p  is  assumed  to  equal  the  radiation  lost  per  unit  volume  from  an 
infinite  cylinder  of  equivalent  T,  p  and  R.  The  reabsorption  effects 
included  are  those  of  the  cylinder  at  uniform  temperature  and  pressure. 
Only  two  extremes  of  reabsorption  effects  are  recognized  in  the  final 
calculations  of  power  lost:  the  radiation  of  a  given  frequency  is  treated 
as  being  either  transparent  or  black. 

Our  calculations  show  that  radiation  losses  rise  rapidly  with 
temperature,  peaking  at  approximately  18,  000  K  where  ionization  is 
complete.  It  is  important  therefore  to  know  the  operating  temperature 
accurately,  A  small  error  in  determining  the  operating  temperature  can 
lead  to  a  large  increase  in  radiation  lost.  Since  the  operating  temperature 
is  intimately  related  to  the  propagation  mechanism,  it  is  also  important 
to  understand  thoroughly  how  LSC  waves  propagate. 

Black  body  surface  radiation  is  less  important  than  transparent 
radiative  losses  for  the  configurations  studied.  However,  the  black  body 
radiation  is  not  negligible,  and  its  importance  increases  with  plasma 
diameter.  Future  work  should  include  modeling  of  the  radial  temperature 
gradients  so  that  better  estimates  of  the  surface  radiation  lost  can  be 
made. 
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6.  ESTIMATE  OF  LASER  SUPPORTED  COMBUSTION  WAVE  PROPERTIES 


To  enable  the  flowing  hydrogen  core  of  the  laser -heated  rocket  to 
absorb  laser  energy,  it  must  be  heated  to  a  temperature  where  thermal 
ionization  takes  place,  i.e.,  to  temperatures  of  order  10,  000  K.  In  steady 
state  the  mechanism  for  providing  this  heating  is  a  laser  supported  combustion 
(LSC)  wave,  which  is  a  zone  of  rapid  temperature  increase  from  ambient 
temperature  to  temperatures  where  absorption  of  the  laser  energy  can  take 
place.  Laser  energy  flows  through  this  zone  from  the  cold  side  and  is 
absorbed  on  the  hot  side.  The  hot  plasma  transmits  energy  forward  to  the 
cold  gas  by  conduction  and  by  radiation  emission  from  the  plasma,  this 
providing  the  energy  necessary  to  heat  the  cold  gas  up  to  temperatures  where 
it  can  absorb  laser  energy.  Thus  the  wave  is  supported  by  the  laser  energy, 
which  indirectly  converts  cold  flowing  gas  to  hot  flowing  gas,  using  laser 
absorption  in  the  hot  gas,  plasma  emission,  and  conduction  and  plasma 
radiation  absorption  in  the  cold  gas  as  mechanisms. 

Such  waves  have  been  observed  in  air,  in  a  non-steady  geometry, 
by  focusing  laser  intensities  of  10^  to  10^  W/cm^  on  solid  targets. 

Evaporation  of  the  targets  provides  electrons  which  initiate  absorption  of 
laser  energy  in  the  gas  in  front  of  the  target,  and  a  temperature  wave  is 
observed  to  move  away  from  the  target  toward  the  laser  beam,  at  subsonic 
speeds.  This  wave  is  an  LSC  wave  moving  relative  to  the  still  air.  It  is 
observed  to  move  at  a  definite  speed,  which  depends  on  the  laser  intensity, 

(At  intensities  greater  than  10^  W/cm^,  the  rapid  energy  deposition  in¬ 
creases  the  pressure  greatly  and  caused  shock  waves  to  precede  the  heated 
zone,  a  phenomenon  called  a  laser -supported  detonation  (LSD)  wave  ), 

The  laser-heated  rocket  flow  is  conceived  to  begin  with  a 
stationary  LSC  wave  providing  the  mechanism  to  heat  the  flowing  hydrogen. 

The  wave  is  made  stationary  by  flowing  parallel  to  the  laser  beam  at  just 
the  speed  with  which  the  wave  would  want  to  propagate  towards  the  beam. 

This  wave  forms  the  initial  condition  for  the  flow  in  the  plasma  core,  so 
it  is  necessary  to  know  its  speed  for  a  given  laser  intensity,  to  establish 
the  rate  of  mass  flow  into  the  core. 


Experimentally,  LSC  waves  in  air  have  been  observed  by  the  non¬ 
steady  method  described  above  and  curves  of  intensity  vs.  speed  derived 


(Ref.  28).  Theoretically, 

•  T  29,  30,3r 

air  LSC  waves  ’  , 


a  number  of  attempts  have  been  made  to  model 
although  their  success  in  comparing  with  experi¬ 


ment  has  been  limited.  There  is  no  information,  either  experimental  or 


theoretical,  on  LSC  waves  in  hydrogen.  Since  radiation  from  the  gas  plays 


a  crucial  role  in  LSC  wave  properties,  and  since  radiation  properties  of 
hydrogen  are  quite  different  from  those  of  air,  one  cannot  easily  extrapolate 
from  air  to  hydrogen.  The  lack  of  reliable  information  on  the  intensity  vs. 
speed  relation  of  LSC  waves  in  hydrogen  must  be  considered  one  of  the 
major  uncertainties  in  the  study  of  the  flowing  core  laser-heated  rocket. 

Nevertheless,  it  was  necessary  to  estimate  the  intensity  vs. 
speed  curve  for  hydrogen  in  order  to  provide  initial  conditions  for  the 
plasma  core  streamtube  calculation.  For  this  purpose,  the  simple  model 
of  Jackson  and  Nielsen^^  was  adapted  to  hydrogen.  In  this  model  we  needed 
to  characterize  the  high  temperature  emission  and  low  temperature  absorp¬ 
tion  of  hydrogen.  For  the  latter,  a  classical  continuum  absorption  coefficient 
was  assumed.  For  the  former,  the  radiated  power  was  taken  to  be  the 
continuum  radiated  power  per  unit  volume,  as  found  by  Yos  in  Ref.  32. 

With  these  assumptions  the  Jackson-Nielsen  model  was  used  to  obtain  the 
intensity  vs.  speed  curve  for  hydrogen  at  1  atm  and  3  atm. 

The  results  are  shown  as  the  solid  lines  on  Fig.  6.  1.  Also  shown 
there  is  the  Jackson-Nielsen  result  for  air  at  1  atm  as  the  dashed  line,  and 
the  data  for  1  atm  air  of  Ref.  28.  It  is  seen  that  the  hydrogen  estimates  do 
not  differ  much  from  the  air  theory  in  the  range  4E5  to  1E6  W/cm^.  Also 
evident  is  the  disagreement  between  the  air  theory  and  experiment,  with 
the  theory  about  a  factor  of  3  lower  than  the  experimental  results. 

In  spite  of  the  uncertainty,  a  point  near  the  3  atm  line  at  3.  67E5 
W/cm^  (3.67E9  W/m^)  and  1340  cm/s  (13.4  m/s)  was  used  for  the  initial 
conditions  of  the  core  calculation.  This  corresponds  fo^  T^  =  300  K  to  a 
mass  flow  per  unit  area  of  0.3316  g/cm^-s  (3,  316  kg /m  -s). 

For  the  10  atm  core  calculations,  it  was  assumed  that  the  inten¬ 
sity  vs.  mass  flow  per  unit  area  was  the  same  as  for  3  atm;  this  reduces 
the  speed  to  400  cm/ s  (4.  0  m/ s). 


These  estimates  were  made  before  the  work  described  in  Sec.  5  was 
undertaken  to  define  hydrogen  radiation  emission.  That  work  could  be 
used  to  develop  a  much  more  sophisticated  model  of  LSC  waves  in  hydrogen, 
which  avoids  most  of  the  assumptions  inherent  in  the  model  of  Jackson  and 
Nielsen,  and  which  uses  a  much  more  accurate  description  of  radiation 
emission  and  absorption  than  is  possible  with  the  formulas  of  Ref.  32. 
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7.  FLOWING  CORE  STREAMTUBE  MODEL 


7.1  INTRODUCTION 

The  central  streamtube  of  the  flowing-core  laser-heated  rocket 
combines  a  number  of  physical  phenomena.  It  absorbs  laser  radiation, 
it  dissociates  and  ionizes,  it  emits  plasma  radiation,  and  it  flows  through 
a  variable  area  channel.  All  these  phenomena  must  be  included  in  a  model 
to  calculate  the  properties  of  the  streamtube.  The  analysis  described  in 
the  final  report  on  the  previous  contract  (Ref.  1)  included  laser  radiation 
absorption  in  a  variable  area  channel,  but  treated  the  gas  as  perfect,  with 
constant  specific  heats,  and  did  not  include  loss  of  energy  by  plasma  radia¬ 
tion.  Hydrogen  is  fully  dissociated  at  about  8000K  (with  dissociation  energy 
2.  16  X  10  J/kg)  and  is  50%  ionized  between  16,  000  and  17,  OOOK  (with 
ionization  energy  13.  1  x  10®  j/kg),  while  the  total  enthalpy  change  from 
300K  to  8,  OOOK  is  3.  76  x  10  J/kg  and  from  8,  000  to  16,  OOOK  is  8.  37  x  10^ 
j/kg.  Thus  it  is  clear  that  the  dissociation  and  ionization  processes  are 
important  energy  sinks,  and  hydrogen  should  be  treated  as  a  real  gas.  The 
loss  of  energy  by  plasma  radiation  emission  both  reduces  the  energy  in  the 
streamtube  and  defines  an  energy  loss  against  which  the  walls  must  be 
protected,  so  it  is  also  important  to  include  this  physical  phenomena. 

A  model  of  the  flowing-core  streamtube  was  constructed  which 
included  variable  area,  one-dimensional  flow  of  dissociating  and  ionizing 
hydrogen  in  chemical  equilibrium,  with  laser  energy  absorption  by  inverse 
Bremsstrahlung,  and  plasma  radiation  emission.  This  model  was  reduced 
to  the  solution  of  three  first  order  differential  equations  which  were  solved 
by  straightforward  integration  on  a  computer.  This  section  will  describe 
the  model,  and  give  the  results  of  some  calculations  for  configurations  of 
interest. 

7.2  BASIC  EQUATIONS 

The  one-dimensional  mass  conservation  equation  for  a  gas  of 
density  p,  speed  u  in  a  channel  of  area  A  gives  the  mass  flow  rate  as 


rh 


(7.1) 
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The  momentum  conservation  equation  relates  the  pressure  gradient  to  the 
velocity  gradient  by 


^  +  DU  =0  (7.2) 

dx  dx 

The  energy  equation  has  laser  absorption  terms  and  plasma  emission  terms. 
The  conservation  of  energy  requires  the  divergence  of  the  total  energy  flux 
to  vanish.  This  total  energy  flxix  has  three  parts  in  the  present  problem. 

One  is  the  total  enthalpy  flux  of  the  gas,  p  uAH,  where  H  is  the  total  enthalpy. 
The  second  is  the  laser  energy  flxox,  AI,  where  I  is  the  laser  intensity  per 
unit  area,  and  the  third  is  the  plasma  emission  radiant  energy  flux  vector 
S.  The  energy  equation  is  then 

R 

(p  u  AH  +  AI)  +  J"  2tt  V  V  *  S  dr  =  0  (7.  3) 

0 

The  laser  energy  is  absorbed  with  absorption  coefficient  accord¬ 
ing  to 


d  AI 
dx 


k  AI 
Xj 


whose  solution  is,  with  (  denoting  conditions  at  x  =  0, 


(7.4) 


AI  =  (AI)j  e 


T  = 


/ 


k  dx 

Lj 


(7.5) 


Use  of  Eqs.  (7.  4),  (7.  5)  and  (7.  1)  in  (7.  3)  allows  the  energy  equation  to  be 
written  a  s 


+  u^  ^  ^ 
dx  dx  dx  dx 


dQ 

dx 


(AI).  e 


R 


2tt  r  V  * 


0 


S  dr 


ih 


(7.6) 


(7.  7) 
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This  is  written  so  the  laser  absorption  term  is  recognized  as  a  heat  addition, 
and  the  plasma  emission  term  as  a  heat  loss. 

The  problem  is  then  one  of  heat  addition  to  a  one-dimensional  nozzle 
flow,  with  the  heat  addition  term  coupled  to  the  flow,  since  kj^  and  V  •  S 
depend  on  the  local  properties  of  the  gas. 

The  conventional  way  to  solve  this  problem,  as  described  in  Ref.  1, 
is  to  take  the  area  A(x)  as  given,  and  solve  for  dp,  du  and  dp  from  (7.  1), 

(7.  2)  and  (7.  6).  The  results  show  the  sonic  singularity  at  M  =  1,  which 
makes  for  difficulties  in  the  numerical  integration  of  the  equations  when 
there  is  heat  addition  or,  real  gas  effects,  both  of  which  are  present  here. 

An  alternative  method  of  solution  was  therefore  used,  which  led  to  simple 
numerical  solutions. 

Instead  of  taking  the  area  as  given,  we  take  the  velocity  distribution 
u(x)  as  known.  This  leaves  the  thermodynamic  variables  p,  p,  T  to  be 
found  from  (7.  2)  and  (7.  6)  and  the  equation  of  state.  Then  A(x)  is  found 
from  (7.  1),  since  the  mass  flow  rate  m  is  given.  This  is  a  perfectly 
acceptable  method  of  solving  the  problem,  as  long  as  one  chooses  u(x)  so 
the  resulting  area  distributions  are  realistic,  and  this  can  be  accomplished 
■with  a  little  experience. 

Equations  (7.  2)  and  (7.  6)  can  be  written  easily  in  terms  of  dp  and 
dT  by  using  the  equation  for  the  enthalpy  of  equilibrium  hydrogen  in  the 
form 


h  =  h  (p,  T) 


=  (It)  •  (It)„  = 


(7.8) 


If  this  is  inserted  into  (7.  6),  and  (7.2)  used  to  eliminate  dp,  (7.  2) 
and  (7.6)  become 


dp  du 

dx  “  "  ^  ^  dx 


(7.  9a) 


dT 

dx 


(7.  9b) 
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when  the  thermodynamic  properties  of  hydrogen  and  the  radiation  properties 
which  are  contained  in  dQ  are  defined,  the  specification  of  u(x)  enables  a 
straightforward  integration  of  these  two  eqiiations  from  any  given  initial 
condition.  From  their  solution  and  the  gas  law,  the  radius  of  the  stream- 
tube  is  obtained  from  (7. 1)  as 

R  =  (m/TTpu)^/^  (7.10) 

Actually  the  emisjsion  term  in  (7*  7)  involves  R,  so  (7*  10)  is  used  for  R 
during  the  course  of  the  integration# 

7.  3  THERMODYNAMIC  PROPERTIES  OF  EQUILIBRIUM  HYDROGEN 


We  need  the  derivatives  of  enthalpy  with  respect  to  p  and  T,  the 
equation  of  state  which  relates  p  and  T,  and  the  composition  of  the  mixture 
of  hydrogen  molecules,  atoms,  ions  and  electrons.  The  pressure  deriva¬ 
tive  of  h  can  be  expressed  in  terms  of  the  equation  of  state  by  a  classical 
use  of  reciprocity  relations  between  the  derivatives  of  h  and  entropy,  which 
are  in  turn  obtained  from  the  second  law#  These  lead  to 


so  that 


(7.11) 


We  characterize  the  hydrogen  as  a  mixture  of  perfect  gases  whose 
components  have  number  densities:  molecules  atoms  n^,  ions  nj  and 

electrons  n^.  We  do  not  expect  to  get  to  a  temperature  where  other  species 
are  significant,  nor  are  there  multiply  ionized  atoms,  so  nj  =  n^#  As 
reaction  coordinates  we  use  the  fraction  of  molecules  dissociated  0  and  the 
fraction  of  atoms  ionized,  a  : 


(Ha  +  nj)/2 

"m  <"a  ^ 


a 


In  terms  of  the  original  number  of  particles  per  unit  volume. 


(7.12) 


n 

o 


=  n 


M 


(^A  +  nj)/2 


the  number  density  of  the  species  are 


(7.  13) 


n  =  n  (1  -  B)  , 
M  o 


Ua  =  2n^  B  (1  -  a)  .  ^^i  =  ^E  =  2n^Ba  (7.14) 
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and  the  corresponding  partial  pressure  are  found  from 


p.  =  n.  kT  (7.15) 

where  k  is  the  Boltzmann  constant.  The  partial  mass  densities  are,  in 
terms  of  the  particle  masses  m^, 

and  the  total  number  and  mass  densities  are 


n 


n  [1  +  B  (1  +  2a)]  ,  p  =  n  m, , 
o  o  M 


(7.  17) 


Whenever  convenient  we  will  use  the  facts  that  to  an  excellent  approximation 
for  hydrogen,  m^  =  mj  =  mj^/2. 

The  equation  of  state  is  the  sum  of  the  partial  pressures. 


p  =  pZRj^T.  Z  =  l+p{l  +  2a),  Rj^  =  k/mj^  (7.18) 

where  Z  is  the  compressibility  factor  and  Rj^  the  gas  constant  for  the  mole¬ 
cules.  We  have  now  expressed  p  in  terms  of  p,  T  through  a  and  8  by  (7.  18). 
Since  p  and  T  are  the  dependent  variables  in  the  basic  differential  equations 
(7.  9),  we  need  a,  8  in  terms  of  T,  p. 

Since  we  deal  with  hydrogen  in  thermochemical  equilibrium,  a  and 
8  are  determined  by  the  Law  of  Mass  Action  in  terms  of  the  thermodynamic 
variables  of  the  mixture.  There  are  two  reactions  occuring,  dissociation 
and  ionization: 


^  2H  H  ==  +  E‘ 

For  each  the  equilibrium  constant  provides  a  relation  between  the 
partial  pressures  and  the  partition  functions  of  the  species,  involving  the 
heat  of  reaction*  These  relations  are 
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^1 


2TT  m. 


_,5/2  ^ell^elE  -h_°  m,/kT 
(kT)  - - -  e  D  I 

^elA 


(7.  19b) 


where  h  is  Planck's  constant,  and  are  the  vibrational  and  rotational 
partitiolP functions  of  Qeii  are  the  electronic  partition  functions  of  the 
species,  hj^°  is  the  heat  of  dissociation  of  a  molecule  per  unit  mass  and 
hi°  is  the  heat  of  ionization  of  an  atom  per  unit  mass.  From  (7.  14),  (7.  15), 
(7.  17)  and  (7.  18)  the  partial  pressure  ratios  are 


Za  _  4B^  (1  -  g)^ 


n  kT 
o 


2  2 
=  4B  (1  -  g)  JS. 

1  -  B  Z 


(7.20a) 


P  P  .  i->2  2 

^E  _  4B  g 

Pa  2  B  (1  -  0.) 


n  kT 
o 


oo  2 

2  B  a  p 
'l  -  a  Z 


{7.20b) 


The  rotational  partition  function  may  be  approximated  by  T/2  0  where 
0  =  87.  62K  for  H2  and  the  factor  of  2  comes  from  the  symmetry  of  the 

niolecule.  This  is  a  valid  approximation  for  T  »  0^;  it  is  about  10%  low 
at  300K  but  only  0.  5%  low  at  lOOOK.  The  vibrational  partition  function  is 

Q  =  (1  -  e  where  0^  =  5983K  for  H2.  The  electronic  partition 

function  for  the  rriolecule  may  be  replaced  by  its  ground- state  statistical 
weight  1,  since  its  first  excited  state  lies  very  high  (132,  OOOK).  For  the 
electron  QelE  “  ^  because  of  the  two  electron  spin  states,  and  for  the  ion, 
there  are  no  internal  states,  Qgjj  =  !• 
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For  the  atom 


^elA 


.2 

3 


exp  (h^  m^/kT  j  ) 


(7.21) 


The  first  term  is  the  usual  grbund- state  weight  2,  The  subsequent  terms 
are  needed  here,  even  though  hj  m^/k  =  158,  OOOK,  because  we  may  be 
interested  in  hydrogen  atoms  up  to  20,  OOOK.  The  upper  limit  cut-off  of 
the  sum  can  be  obtained  in  various  ways,  and  may  depend  on  temperature. 
Here  we  took  it  constant  at  17,  since  the  effect  of  the  sum  is  rather  small. 
If  (7.20)  and  (7.21)  are  used  in  (7.19)  we  have 


2  2 
B  (1  -  g) 

(1  -  B)  Z 


t3/2 

—  -  V 


(7.22a) 


2  B 

(1  -  a)  Z 


2Tt  k 

.  I 


(7.22b) 


m^/k 


These  two  relations  express  a,  P  in  terms  of  p,  T,  as  desired,  and  could  be 
used.  However,  they  are  algebraically  complicated,  and  difficult  to  solve 
for  a,  P.  For  hydrogen  they  can  be  considerable  simplified  with  little  loss 
inaccuracy.  Hydrogen  in  equilibrium  dissociates  almost  completely  before 
it  ionizes  significantly.  We  can  thus  separate  the  two  reactions  and  take 
a  =  0  in  the  first  one  while  P  goes  from  0  to  1.  Then  in  the  second,  we 
can  take  P  =  1  while  a  increases  from  zero.  These  two  are  separated  by 
a  temperature  T  below  which  there  are  only  molecules  and  atoms,  and 
above  which  there  are  only  atoms,  ions  and  electrons.  Additionally,  in  the 
non-ionizing  region,  the  temperature  is  too  low  for  the  higher  electronic 
states  of  atom  to  contribute,  so  the  squared  term  on  the  right  of  (7.22a) 
can  be  taken  as  unity.  Then  we  have 
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,  a  =  0 


(7.  23a) 


«  'd<^> 

T  <  T  :  — — T  =  - 

1  P 


(Vim  ^  k 

‘'p 


3/2 


2kT  T 
r 


2  f,  (T 

*  CC  I  n  ^ 

T  >  T  ;  - r  =  -  ,  B  =  1 

X  1  2  p 

1  -  a 


(7.23b) 


‘l'^>  -[  —  2 


kT^/2 


1 


2  2 

These  are  both  linear  in  B  or  a  ,  and  yield 


T<T*:  a  =  0,  Z  =  l+B,  B=(l+  p/^q) 


-1/2 


(7.  24a) 


T  >  T 


B  =  1,  Z  =  2  (1  +a), 


a  =  (1  +  p/fj)‘^^^ 


(7.  24b) 


which  are  very  simple  expressions  for  oc,  P  in  terms  of  p  and  T,  the  pri¬ 
mary  dependent  variables. 

From  (11),  we  need  the  derivative  of  In  p  with  In  T  at  constant 
pressure,  which  from  the  gas  law  (7.  18)  is 


h  In  p  \ 

a  In  T  / 
'V 


1 


(  ^  z\ 


(7.25) 


The  derivate  of  Z  is  easily  foiand  from  (7.  24)  in  terms  of  the  derivatives 
of  a,  B,  which  are  obtained  by  differentiating  (7.  24)  and  using  the  definitions 
of  fjj,  fj  given  in  (7.23).  The  results  are 
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(7.  26a) 


(7.  26b) 


The  remaining  thermodynamic  property  needed  is  Cp,  the  tempera¬ 
ture  derivative  of  the  enthalpy.  The  translation  specific  heats  of  the  atoms, 
ions  and  electrons  are  5k/2  per  particle.  This  is  the  total  contribution  for 
the  ions  and  electrons,  which  thus  have  a  specific  heat  per  unit  mass  of 


c  _  =  5k/2m  ,  c  =5k/2m  =c  m/m  (7.27) 

pi  I  pE  E  pi  I  E 

Their  enthapies  are  the  integrals  with  respect  to  T,  plus  the  dissociation 
energy  hjj°  and  the  ionization  energy  hi°  for  the  ions 

hf  =  Cpj  T  +  hj^  +  hj  ,  h^  =  Cp^  T  (7.  28) 
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The  atoms  have  a  translational  contribution  5kT/2m^,  carry  the  dissocia¬ 
tion  energy  hp°,  and  also  have  a  contribution  from  the  electronic  states 
derived  from  Qg^A. 


,  _2  S  In  Q  ,  ^ 

.  _  k  T  _ elA  o  I  _2 

^elA  m^  ST  ^  \  ’  ^1 


where  83  has  been  defined  in  (7.  26b).  The  resulting  atom  enthalpy  is 


.  SkT^.o^.o/  ®2 

h  -  - -  +  h  +  h  I  1  -  „ 

A  2  m^  D  I  y  Sj^ 


(7.29) 


Finally  the  molecules  have  a  translational  contribution  5kT/2mj^,  a 
rotational  contribution  2kT/2mj^  for  the  two  rotational  degrees  of  freedom, 
corresponding  to  =  T/2  9^,  and  a  vibrational  contribution  found  from 

-0  /T  "V 

Q  =  1  -  e  as 

V 


h 

V 


kT 


m 


M 


where  is  defined  in  (7.26a).  Thus  the  molecular  enthalpy  is 


h 

m 


7kT 


2  m 


M 


kT 

m 


M 


(7.30) 


The  enthalpy  of  the  mixture  is  the  sum  of  the  species  enthalpies 
weighted  by  their  mass  fractions  in  the  mixture: 


h  =E  h.  p./p  (7.31) 

The  mass  fractions  are  found  in  terms  of  a  and  9  by  (7.  14),  (7.  16)  and 
(7.17): 

P  Pa  Pj  Pe 

=  -A=B,i.a),  — =ea  —  (7.32) 
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Using  (7.  28),  (7.  29),  (7.  30)  and  (7.  32)  in  the  enthalpy  expression  (7.  31), 
we  find  the  enthalpy  for  this  mixture  as 


h  =  c  ^  T  I  Ml  +a)  +  (1  -  B)  (0.7  +  0.2  ilr  ) 

pi  (  V 

+  0.4B  [Sjj/ZT  +  (a+ (1  -  a)  t^j^)  Gj/T  j| 


(7.  33) 


c  .=  5k/2mj.  ’^elA  =  ^  -  <^1/^2) 


The  equilibrium  c  is  obtained  by  differentiating  h  with  respect  to 
P 


T  at  constant  p; 


c  =  c  _ 

p  pi 


h  < 

(1  -  B)  (0.7  +  0.2cp  )  +  B  (1  +  a)  +  0.4  B  (1  -  a)  —  —  ^  — 

\^l 


+  c  ^  T 

pi 


\aT/ 


u  w 

0.3  -  0.2  +a  +  0.2-^+  0.4^(a  +  (1  -  a)  ^ 


+  c  TB 
pi 


I  — \ 

Ut/ 
p  L 


(7.34) 


cp  =  (6  /T)  e 

V  V 


j=i 


Equations  (7.  11),  (7.  18),  (7.23),  (7.24),  (7.25),  (7.  26)  and  (7.34) 
express  p,  (bh/bp),^,  and  Cp  in  terms  of  p,  T,  and  thus  express  all  the 
unknowns  except  dQ  in  the  basic  differential  equations  (7.  9)  in  a  form 
usable  for  numerical  integration 
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A  check  on  the  thermodynamic  properties  used  can  be  obtained  by 
comparing  some  calculations  made  using  the  present  formulas  with  those 
of  Patch  (Ref.  33).  (The  enthalpies  of  Ref.  33  have  had  hjj°  =  2.  16  E  5  j/g 
added,  since  in  that  reference  the  H  atom  is  taken  as  the  reference  for 
enthalpy,  while  here  we  have  taken  the  Ho  molecule.  )  The  comparison  is 
made  in  Table  7.  1  for  P,  a  and  h,  with  T*  =  8000K.  Comparison  shows  the 
errors  to  be  mostly  a  few  percent,  compared  to  the  very  elaborate  scheme 
of  Ref.  33.  The  worst  errors  are  at  high  values  of  T  and  p.  The  applica¬ 
tions  in  the  presefat  report  do  not  exceed  16,  OOOK  at  10  atm,  and  are  all 
above  8,  OOOK,  so  the  approximations  appear  satisfactory.  The  major 
source  of  error  is  probably  the  use  of  17  terms  in  the  sums  of  Sj  and  S2, 
instead  of  a  variable  cut-off. 

7.4  ABSORPTION  OF  LASER  ENERGY 

The  first  term  in  the  heat  addition  expression,  Eq.  (7.  7),  gives  the 
absorption  of  the  incident  laser  energy  by  the  plasma.  It  is  dependent  on 
k  ,  the  absorption  coefficient  of  the  plasma.  Based  on  the  studies  of  Ref.  1, 
w^  use  as  the  absorption  mechanism  inverse  Bremsstrahlimg,  both  electron- 
ion  and  electron-neutral.  The  expressions  for  kj_^  have  been  given  in  Ref.  1. 
The  electron-ion  absorption  is  given  in  Eq.  (2.  1)  of  Ref.  1  as 

where  is  the  electron-ion  absorption  cross-section  and  v  is  the  laser 
frequency  (for  10.6p.m  in  the  present  case).  The  cross-section  is  given 
in  Eq.  (2.  2)  of  Ref.  1  as 


El 


4 

3; 


(^) 


1/2 


6 

Z  e 


1.63  X  10 


-32 


h  cm„  v' 

P  ilj 


,1/2 


(cm^) 


where  Ze  is  the  ion  charge,  e  is  the  electron  charge  and  c  is  the  speed  of 
light.  Combining  these  two  equations  we  find 


where  T  is  in  K  and  n^,  n^  are  in  cm 


-3 
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TABLE  7.  1 


COMPARISON  OFa,  g,  h  WITH  REFERENCE  33 


p  (atm) 

T  (K) 

B  (Ref.  2) 

%  Error 

1 

2000 

8.  34E-4 

8.  lOE-4 

3.  0 

10 

2000 

2.64E-4 

2.  56E-4 

3.  1 

1 

4000 

0.  645 

0.622 

3.8 

10 

4000 

■  0.258 

0.  243 

6.0 

1 

6000 

0. 99389 

0.  99251 

0.  14 

10 

6000 

0. 9435 

0. 9319 

1.2 

1 

8000 

0.  99994 

0.99926 

7.3 

10 

8000 

0. 9945 

0.  9927 

1.  9 

p  (atm) 

T  (K) 

g 

a  (Ref.  2) 

%  Error 

1 

10,000 

0.2126E-1 

0.2192E-1 

-3.  0 

10 

10, 000 

0.  672E-4 

0.  706E-4 

-4.  8 

1 

14, 000 

0.2922 

0.  3050 

-4.2  . 

3 

14, 000 

0.  174 

0.  185 

-5.9 

10 

14, 000 

0. 0962 

0.  105 

-8.4 

1 

20, 000 

0.  895 

0.  922 

-2.9 

10 

20,  000 

0.  535 

0.651 

-17.  8 

p  (atm) 

T  (K) 

h  (J/g) 

h  (j/g)  (Ref.  2) 

%  Error 

3 

12,000 

5.  604E5 

5.  57E5 

0.5 

10 

12, 000 

5.  199E5 

5.  16E5 

1.  1 

3 

14, 000 

8.  123E5 

8.  03E5 

0.  8 

10 

14, 000 

6.  908E5 

6.  74E5 

2.  5 
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The  electron  and  ion  number  densities  are  related  to  the  density  by 
(7.  14)  and  (7.  17)  as 


=  nj  =  2  e  a  o/mj^ 


The  above  expression  for  is  a  semi-classical  result.  To 

account  for  quantum  -  mechanical  effects  we  multiply  by  a  Gaunt  factor  of 
1.  5,  which  is  an  average  value  for  the  range  of  interest.  The  final  expres¬ 
sion  for  the  electron-ion  inverse  Bremsstrahlung  absorption  coefficient 
is  then 


‘'lei 


=  8.  7  X  10^ 


o}  3^  1357/T 

-  (e 

JT 


1) 


(7.35a) 


a  =  15  for  cgs  units,  a  =  11  for  mks  units. 

Of  course,  this  expression  is  only  valid  where  there  are  ions,  a  0,  so 

_ 

in  our  approximation  for  the  chemistry,  it  is  to  be  used  only  above  T 
where  3=1. 

The  absorption  coefficient  for  electron-neutral  inverse  Bremsstrah¬ 
lung  is  also  given  in  Ref.  1,  Eq.  (2.  6),  as 


k..  =  Q-r^nvT 

XjE  n  e  n  e  n 


[l-exp  {-hpV 


where  is  the  absorption  cross-section  and  n^g-  is  the  neutral  number 

density.  Equation  (2.  8)  of  Ref.  1  gives 


-45 

2. 96  X  10  T 
1  -  exp  (-hpV/kT) 


where  I  may  be  approximated  by  the  first  term  of  Eq.  (2.  9)  of  Ref.  1, 
when  */T/  158,  000  is  small,  as 


1/2 

5=  4.  862  (kT/mj  h°) 
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Since  for  hydrogen  mj  hj°/k  =  158,  000  K,  then  we  have 


^LEN  ~  1*6x10  ^pF  n^  n^^  exp  (-0.  0122  Vt  )  (cm 
where  n^^,  n^^  are  in  cm”^. 

The  niunber  densities  are  found  from  (14)  and  (17),  where  we 
recognize  that  our  chemical  model  has  8  =  1  whenever  there  are  electrons 
present,  so 


=  "n  °  "a  =  ^  “"'m 

Then  the  final  expression  for  is 


k  =  5.  7  X  10 
LEN 


a  (1  -  a)  /T  exp  (-0.  0122  //t" 


(7.  35b) 


b  =  9  for  cgs  units  b  =  5  for  mks  units 

The  expression  for  kj_^  to  be  used  in  the  first  term  of  (7.  7),  and  in 
the  definition  of  optical  depth  T,  Eq.  (7,  5),  is  the  sum 


^LEI  ^LEN 


(7.36) 


A  plot  of  (dashed  lines)  and  kLj.j^  (solid  lines)  is  given  in  Fig.  7.  1  as 

kij/p  vs  T.  It  can  be  seen  that  electron-ion  absorption  dominates  at  temp¬ 
eratures  above  11,  OOOK  at  10  atm,  and  above  9,  OOOK  at  1  atm.  The  sum 
kj^  follows  the  kj  -p^j  curves  closely  above  12,  OOOK, 

7.  5  PLASMA  RADIATION  EMISSION 


The  second  term  in  dQ,  Eq.  (7.  7),  is  the  loss  of  energy  caused  by 
emission  of  radiation  from  the  plasma 


R 

J  2tt  r  V  ♦ 
0 


S  dr 
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The  model  used  for  this  term  is  discussed  in  the  section  of  this  report 
entitled  Radiative  Emission.  The  model  is  one  of  a  uniform  cylinder  of 
plasma,  so  the  emission  at  any  station  is  expressed  in  terms  of  a  volume 
emission  and  the  cross-sectional  area 

R 

y*  2TT  r  V  '  S  dr  =  F  R^ 

0 

The  general  form  of  P,p,  as  deduced  in  the  section  on  Radiative 
Emission,  from  Eq.  (5.32),  is 


where  is  the  atom  ground  state  population,  in  cm~  ,  and  the  constants 

Cj  and  C2  depend  on  the  pressure  level  and  the  relative  importance  of  line 
and  continuum  radiation.  The  relation  between  and  n^  is  through  the 
electronic  partition  function  of  (7.21),  whose  first  term  relates  to  the 

ground  state. 


n 


AG  A 


Q 


elA 


Sj  exp  (-  hj  m^/kT) 


When  n^  is  expressed  in  terms  of  a,  P  and  p  from  (7.  14)  and  (7.  17),  the 
expression  for  P^  is,  in  mks  units. 


=  5.1 


X 


e  (1  -  g)  p 

^1 


© 


(7.  38) 


The  values  of  the  constants  C^,  depend  on  the  pressure  level  and  laser 
power,  as  described  in  the  section  on  Radiative  Emission  from  Hydrogen. 
The  particular  values  used  for  the  cases  calculated  will  be  given  when  the 
results  of  the  calculations  are  discussed. 


7.  6  SPECIFICATION  OF  VELOCITY  DISTRIBUTION 


As  described  above,  in  the  subsection  on  Basic  Equations,  the 
velocity  distribution  in  the  plasma  core  streamtube  is  chosen  as  the  indepen¬ 
dent  parameter  rather  than  the  area  distribution,  since  it  greatly  simplifies 
the  calculation  of  the  flow,  by  avoiding  the  explicit  appearance  of  the  singu¬ 
larity  at  the  station  at  which  M  =  1. 
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Several  forms  of  u  were  tried,  to  see  if  reasonable  nozzle  shales  were 
attained.  The  simplest  and  most  satisfactory  was  a  logarithmic  distribution. 

The  calculations  were  done  with  the  optical  depth  t  ,  rather  than 
distance  x,  as  the  independent  variable.  This  choice  eliminates  the  explicit 
appearance  of  kL  from  the  differential  equations  when  the  plasma  emission 
term  e  is  not  taken  into  account,  as  can  easily  be  seen  from  Eqs.  (7.5),  (7.7) 
and  7.9).  It  was  used  in  Ref.  1 ,  where  e  was  not  considered.  In  the  present 
case  where  C  was  included,  there  is  no  simplification  afforded  by  using  T, 
since  it  will  require  the  explicit  appearance  of  k^  in  the  plasma  emission 
term,  but  it  was  used  to  provide  some  continuity  with  Ref.  1.  However,  in 
the  future,  it  probably  would  be  better  to  use  x  as  the  independent  variable. 

In  terms  of  T,  the  u  distribution  was  taken  as 


which  goes  from  at  T  =  0  (x  =  0),  the  initial  station,  to  U£  at  the  final 
station  T^.  This  is  a  linear  variation  of  In  u  with  T.  The  nozzle  shapes 
obtained  with  this  can  be  quite  reasonable,  as  will  be  seen  when  the  results 
of  the  calculations  are  presented. 

The  derivative  is  easily  obtained  from 


as 


In  u  =  In  u^  +  (t/t^)  In 

u  In  (u  /u. ) 
du  _  _  f  1 

dT  ”  T 


(7.39b) 


7.  7  SUMMARY  OF  THE  MODEL 

All  the  quantities  necessary  for  the  integration  of  the  differential 
equations  (7.  9)  have  now  been  defined.  In  terms  of  T  as  the  independent 
variable,  (7.  9)  becomes,  using  (7.  5),  (7.  7),  (7.  10),  (7.  11),  (7.  18),  (7.  25) 
and  (7.  37): 


dp  _  ~  pu  du  (7.40a) 

dT  "  R^,  ZT  dT 

M 
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dT 

dT 


TT  e 
k  m 


\ 


(7.40b) 


dx  1 

dT  k_  ’ 

JLi 


(7.40c) 


Z  is  given  in  terms  of  a,  P  by  (7,  24),  and  those  equations  also  give  ct*  3 
in  terms  of  p,  T,  Equation  (7.  34)  gives  c  ,  while  the  derivatives  of  a,  P 
and  Z  are  in  (7.26),  With  u  and  du/dT  from  (7.39),  a  forward  integration 
of  the  three  differential  equations  can  be  performed  from  a  starting  station 
T  =  0  at  which  x  =  0,  T  =  T.,  p  =  p.,  u  =  u^,  to  a  final  station  T^,  where 
u  =  U£,  The  mass  flow  rate  constant  m  is  determined  by  the  given  initial 
conditions  and  the  initial  radius  as 


•  ^ 

m  =  TT  R.  p.  u.  =  A.  p.  u. 

Ill  111 


(7.41) 


which  shows  that  only  is  determined  by  the  solution#  Of  course,  if 

the  laser  power  P  is  given,  is  related  to  the  incident  laser  intensity  by 


I  I  =  P/rr  R.^  (7.42) 

o  1 

The  incident  laser  intensity  is  not  the  same  as  the  initial  intensity  !•, 
because  energy  has  been  added  to  the  incident  gas  in  the  LSC  wave  before 
the  stream  tube  calculation  is  started,  and  we  must  be  careful  to  conserve 
energy.  If  the  calculation  is  begun  at  Tj^,  Pj,  u^,  the  initial  total  enthalpy 
flux  is  m  (h^  +  m^/2),  where  h^  is  obtained  from  (7,33).  If  the  gas  is 
supplied  at  Pq,  Tq,  h^,  u^,  then  its  total  enthalpy  flux  is  rh  (h^  +  Uq'^/2), 
and  the  amount  of  energy  added  up  to  the  initial  state  must  have  come  from 
I^.  Thus 


I.  =  I  -  m  (h.  +  u.^/2  -  h 
1  o  All  o 


u  ^/2) 

o 


(7.43) 
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(7,44) 


where  m. 


is  the  mass  flow  per  unit  area, 


=  m/Aj  =  p.  =  p 


o 


u 

o 


If  the  initial  conditions  are  at  a  temperature  of  order  10,  OOOK,  the  reduction 
in  I  ,  caused  by  the  heating  of  the  incident  gas  by  the  LSC  wave,  is  sub¬ 
stantial. 


The  initial  conditions  for  the  streamtube  are  determined  by  the 
properties  of  the  L.SC  wave  which  pre-heats  the  gas.  As  is  customary  with 
LSC  wave  properties,  we  assume  it  causes  no  pressure  change,  but  only 
temperature  and  density  changes,  so  p^  =  p^.  In  addition  the  LSC  wave 
propagates  against  the  flow  at  a  speed  dependent  on  laser  intensity  and  preS' 
sure  level,  so  I  and  p  determine  If  is  also  given  then  p^  and 

u  =  m^/p^  u^  =  rn^/pj.  Finally,  (7,43)  gives  I.,  and  we  have  all  the 
necessary  initial  conditions.  This  scheme  shows  that  the  properties  of 
the  LSC  wave,  plus  a  choice  of  1^,  p^  and  Tq  are  all  the  conditions  neces¬ 
sary  to  integrate  the  streamtube  equations,  once  u^  and  are  chosen.  The 
integration  then  determines  the  distribution  of  all  gas  properties  along  the 
streamtube  including  the  plasma  emission  e,  and  also  determines  the 
streamtube  shape  R/R.,  The  scale  of  the  streamtube  can  only  be  found 
when  the  laser  power  is  also  given,  whence  R^  follows  from  (7,42), 
However,  it  should  be  noted  that  some  idea  of  the  size  of  the  streamtube, 
is  necessary  to  perform  the  integration,  because  its  size  influences  the 
degree  of  opacity  of  the  plasma,  and  thus  influences  the  choice  of  constants 
in  the  expression  for  e. 

Although  it  is  not  used  in  the  integration,  the  Mach  number  is  an 
interesting- parameter  which  can  be  calculated  when  p  and  T  are  known, 
since  they  determine  the  speed  of  sound  a.  For  a  gas  in  thermochemical 
equilibrium,  the  second  law  of  thermodynamics  and  the  expression  of  h  as 
a  function  of  p  and  p  shows  that 

(bh/3p) 

—  — E  . . 

s  p'^  -  (Sh/Sp) 

P 

Conversion  of  these  derivatives  to  h  as  a  function  of  T  and  p,  and  use  of 
(7.11)  shows  that 
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The  second  derivative  has  already  been  expressed  in  terms  of  Z  by  (7.25). 
From  the  gas  law  (7.  18)  it  is  easy  to  see  that 


Thus  the  Mach  number  squared  is 


-  2  2 
2  u  u 
M  =  —  =  —  ,  „  ^ 

2  T  J  R.  .  Z 
a  I  M 


/  a  In  Z  \  /a  In  Z  \ 

\  a  In  p  /_  c  ^ \a  In  T  / 

\  /T  p  ^  p 


(7.45) 


The  T  derivative  of  Z  is  given  in  (7.  26).  Similar  differentiation  of  Z,  a  and  P 
in  (7.  24)  with  respect  to  p  gives 


T  <  T 


*  /a  In  z \  ^  1  /  ap  \ 

'  \a  In  p  1  +  P  \a  In  .p/, 

*  /a  In  z  \ 

'  \d  In  p  L 


1  +  a  \  a  In  p 


P  (1  -  6^) 
2 

q  (1  -  g^) 
2 


(7.46a) 


(7.46b) 


Thus  the  Mach  number  can  be  calculated  as  the  integration  progresses. 
7,  8  RESULTS  OF  CALCULATIONS 


A  computer  program  has  been  written  to  integrate  Eqs.  (7.40),  and 
has  been  used  for  six  cases  of  varying  laser  power  and  incident  gas  pres¬ 
sure  level.  All  the  cases  were  done  for  the  same  incident  laser  intensity 
of  =  3.67  X  10°  w/m^  (3.67  x  10^  w/cm^)  and  the  same  mass  flow  rate 
per  unit  area,  m^  =  3.  316  kg/m^  (0.  3316  g/cm^),  based  on  an  estimate  of 
the  properties  of  LSC  waves  in  hydrogen,  as  described  in  the  section  on 
LSC  Wave  Properties.  In  all  cases  the  incident  temperature  was  T^  =  300K 
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and  the  initial  temperature  was  T.  =  12,  OOOK.  (It  is  believed  that  the  effect 
of  lowering  by  a  few  thousand  degrees  will  not  change  the  major  results 
of  the  calculations.) 

The  parameters  varied  were  the  laser  power  P  and  the  incident  gas 
pressure  level  Pq  =  p^.  Therefore  the  initial  radius  R.  (and  area  A.)>  and 
the  total  mass  flow  rate  is  also  varied.  The  calculation  matrix  is  given  in 
Table  7.  2. 

The  range  of  power  varied  from  the  experimental  device  of  10  kw  to 
a  5000  kw  device  of  operational  size,  with  an  intermediate  size  of  100  kw. 
The  pressure  of  3  atm  corresponds  to  the  experimental  device,  while  10 
atm  was  chosen  to  show  the  effect  of  operating  pressure.  The  initial 
densities  and  velocities  are  determined  by  p^,  T^  and  as  explained 

above.  In  all  cases,  the  integration  was  carried  out  to  a  final  value  of 
optical  depth  =  4,  with  a  velocity  there  of  u^  =  1.  47  E4  m/s.  The  attempt 
was  to  perform  the  calculations  to  the  throat  of  the  streamtube,  and  in  all 
cases  the  minimum  area  was  in  fact  reached.  The  choice  of  =  4  means 
that  the  laser  energy  s^ill  unabsorbed  at  that  point  is  down  from  the  initial 
value  by  a  factor  of  e  =1,8  E-2,  so  only  2%  of  the  initial  laser  energy  is 
not  absorbed.  From  fhere  one,  there  is  no  longer  significant  heat  addition 
to  the  gas.  There  are  still  losses  due  to  radiation. 

The  constants  C3  and  C2  used  in  Eq,  (7.38)  to  express  the  plasma 
emission  e  were  taken  from  the  material  described  in  the  section  on 
Radiative  Emission,  and  were  appropriate  to  the  sizes  of  plasma  core  for 
each  power  level.  They  are  given  in  Table  7.  3. 

The  properties  of  hydrogen  needed  were  taken  from  standard 
references,  and  are: 


k  =  1.  305  E-23  J/K,  m^  =  m^  =  rn^/2  =  1.  673  E.27  kg 


m^  =  9. 109  E-31  kg,  hp  =  6.  6256  E-34  Js 

6-  =  87.  62  K,  e  =  5983  K 

^  V 


Op  =  52,  OOOK,  Oj  =  158,000  K  , 

The  value  used  for  9^  includes  the  correction  for  anharmonicity. 
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TABLE  7.2 


CASES  CALCULATED 


Case 

P  (kW) 

p  (atm) 
o 

R.  (cm) 

1 

A.  (cm^) 

1 

•  2  2 
m(g/cm  -s  =  10  kg/m  -s) 

1 

10 

3 

9.31  E-2 

2.  72  E-2 

9.03  E-3 

2 

5000 

3 

2.08 

1.36  El 

4.  51 

3 

100 

3 

2.94  E-1 

2. 72  E-1 

9.  02  E-2 

4 

10 

10 

9.31  E-2 

2.72  E-2 

9.  03  E-3 

5 

100 

10 

2.94  E-1 

2.  72  E-1 

9.  02  E-2 

6 

5000 

10 

2.08 

1.36  El 

4.  51 

Case 

TABLE  7.  3 

CONSTANTS  IN  PLASMA  EMISSION 

P(kW)  C3(WM^ 

<=2 

1 

10 

3 

2 

2.4 

2 

5000 

3 

2 

0.4 

3 

100 

3 

1 

2.4 

4 

10 

10 

4/3 

2.4 

5 

100 

10 

3/4 

2.4 

6 

5000 

10 

2.2 

0.4 

-77 


The  integration  of  the  differential  equations  was  performed  at  steps 
in  T  of  0.  05,  so  there  were  80  steps  from  T  =  0  to  4, 

The  profiles  of  T  and  R/R^  vs  x  are  given  in  Figs,  7.2  and  7.  3. 
Figures  7.2a  and  7.2b  are  the  temperature  profiles  for  3  and  10  atm  respec¬ 
tively,  The  location  of  the  throat,  and  the  station  at  which  90%  of  the  initial 
laser  energy  is  absorbed  are  marked.  They  all  show  an  increase  to  a  rather 
flat  maximum,  followed  by  a  drop  to  the  throat,  although  the  throat  is  still 
at  temperatures  of  13,  000  to  15,  000  K.  The  maxima  are  quite  flat,  and 
occur  near  the  90%  absorbed  station.  The  larger  the  power,  the  larger  is 
the  temperature  at  any  station,  as  might  be  expected.  Also,  the  large  the 
power,  the  shorter  the  distance  to  the  throat,  though  this  effect  is  not  great. 
These  distances  vary  from  3.  2  to  4.  0  cm  for  3  atm  and  from  0.  49  to  0.  54 
cm  for  10  atm.  This  points  up  the  major  effect  of  pressure  level.  Because 
the  absorption  coefficient  increases  with  pressure,  the  laser  energy  is 
absorbed  in  shorter  distances  at  the  higher  pressure,  and  the  converging 
part  of  the  streamtube  is  shorter.  The  10  atmosphere  streamtubes  are 
nearly  one-tenth  the  length  of  the  3  atmosphere  streamtubes. 

When  one  combines  the  initial  radii  given  in  Table  7.2  with  the 
throat  lengths  given  in  Fig.  7,2,  the  small  size  of  these  devices  becomes 
evident.  The  lengths  are  determined  by  the  pressure  level  through  the 
absorption  coefficient,  as  just  mentioned.  The  radius  is  determined  by  the 
power  and  the  incident  intensity  I  .  If  1q  of  order  10^  W/cm  is  required  to 
support  the  LSC  wave,  then  radii  of  order  0.  1  to  2  cm  will  follow.  To 
enlarge  the  size,  it  is  crucial  to  operate  at  as  small  values  of  as  possible 
which  means  as  near  as  possible  to  the  threshold  intensity  of  the  LSC  wave. 
This  emphasizes  the  importance  of  having  a  reliable  theory  to  estimate 
the  threshold  of  LSC  waves  in  hydrogen. 

Figures  7.2a  and  7,3b  show  the  normalized  radius  distribution 
R/Rj^  for  3  and  10  atm  respectively.  It  can  be  seen  that  only  small  varia¬ 
tions  of  shape  with  power  occur,  and  the  shapes  seem  very  reasonable. 

Notice  the  vertical  scale  is  not  physical  size,  though  the  horizontal  scale 
is,  so  angles  on  the  figure  are  not  necessarily  correct.  For  the  5  MW 
power  R£  =  2.  08  cm,  so  the  streamtube  shown  on  the  3  atm  figure.  Fig.  7.  3a, 
has  very  nearly  the  same  scales  on  both  axes.  The  same  is  approximately 
true  of  the  100  kW  power  on  the  10  atm  figure.  Fig.  7,  3b,  where  R^^  =  0.  294 
cm.  Other  streamtube  shapes  can  be  produced  by  varying  the  assumed 
velocity  distribution,  and  more  experience  with  calculations  would  show 
what  range  of  shapes  can  be  achieved. 

Notice  that  at  10  atm  the  higher  power  cases  lead  to  fat,  stubby 
streamtubes,  such  as  one  which  is  2,  08  cm  in  initial  radius  and  only  0.  5  cm 
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Fig.  7.2a  Laser-Heated  Streamtube  Temperature  Distribution 
for  3  atm  Initial  Pressure 


*12,000  K 
*3.316  kg/m 
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Laser-Heated  Streamtube  Radius  Distribution 
for  3  atm  Initial  Pressure 


Laser-Heated  Streamtube  Radius  Distribution 
for  10  atm  Initial  Pressure 


long  for  5  MW.  This  may  be  undesirable  geometry  from  flow  criteria, 
and  may  indicate  that  high  pressure  operation  at  large  powers  is  not  a 
useful  combination. 

It  can  be  concluded  that  this  rather  small  program  can  calculate  the 
core  streamtube  quite  easily. 

Other  results  of  the  calculations  are  described  in  the  sections  on 
Radiative  Heating  and  on  Performance. 
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8.  HEAT  LOSSES 


8.  1  INTRODUCTION 

The  flowing  plasma  core  loses  heat  by  emission  of  plasma  radiation, 
and  by  convenction  of  heat  to  whatever  surface  is  in  contact  with  the  hot 
plasma.  The  walls  of  the  rocket  must  be  protected  from  this  heat  load.  In 
order  to  define  what  method  of  protection  should  be  used,  .the  size  of  the 
heat  load  must  be  known.  This  section  is  devoted  to  calculation  of  the  heat 
load  from  both  emission  and  convection. 

The  emission  has  already  beenused  in  the  plasma  core  streamtube 
calculation,  where  it  was  coupled  to  the  plasma  core  temperature  and 
geometry  to  provide  the  energy  loss  term.  The  explicit  calculation  of  the 
amount  of  energy  lost  and  the  resulting  heat  load  is  easily  obtained  from 
the  results  of  the  core  streamtube  calculation,  as  described  in  Subsection 
8.2. 


No  convective  heat  loss  was  taken  into  account  in  the  core  calculation. 
In  order  to  decide  if  this  loss  is  of  sufficient  size  to  require  it  also  to  be 
coupled  to  the  core  calculation,  an  uncoupled  calculation  was  made  by  using 
the  temperature  and  pressure  distributions  found  from  the  core  model. 

These  conditions  were  used  as  external  conditions  in  a  hydrogen  boundary 
layer  calculation  to  estimate  the  heat  loss  to  a  cool  surface  placed  in  contact 
with  the  plasma  core.  If  the  resulting  energy  loss  is  large  compared  to  the 
core  energy,  it  indicates  a  coupling  of  the  convective  loss  and  the  core  flow 
is  needed.  If  the  loss  is  small  compared  to  the  core  energy,  such  a  coupled 
model  is  not  necessary.  In  addition  to  an  impervious  wall  boundary  layer, 
some  cases  of  gas  injection  into  the  boundary  layer  were  calculated  to  see 
if  that  was  a  feasible  method  of  reducing  the  convective  heat  loss.  These 
hydrogen  boundary  layer  calculations  are  described  in  Subsection  8.3. 

Finally,  the  possibility  of  protecting  the  walls  from  the  radiative 
eitiission  by  particles  in  a  buffer  gas  region  surrounding  the  plasma  is 
considered  in  Subsection  8.4. 

8.2  RADIATIVE  HEAT  LOSS 

One  of  the  results  of  the  core  streamtube  calculation  is  the  radiative 
energy  emission  from  the  plasma.  At  any  station,  the  plasma  emission 
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termisTTR^  p^,  as  shown  in  Eq.  (7.37).  The  streamtube  prc^gram 
calculates  both  R  and  P-j  as  function  of  x,  and  the  quantity  TTR  is 

shown  in  Fig.  8.  la  for  3  atm  and  in  Fig.  8.  lb  for  10  atm.  Notice  that  in 
Fig.  8.  la,  the  5  MW  curve  is  plotted  as  TT  R^  P,j./100,  while  in  Fig.  8.  lb 
the  100  kW  curve  is  TTR^  Pj/10  and  the  5  MW  curve  is  tTR  P,j./100. 

These  curves  show  that  the  station  emission  has  a  maximum  some¬ 
what  earlier  than  the  temperature  profiles  shown  in  Fig.  7.  2,  because  R 
decreases  steadily.  At  the  throat,  the  values  are  below  the  initial  values 
because  of  the  small  radius  there. 

In  order  to  estimate  the  heating  received  by  the  wall  or  the  buffer 
gas  by  plasma  emission  we  divide  the  station  emission  by  the  circumference 
2  TTR  of  the  steamtube  at  the  same  station,  to  get  the  distributions  of 

q  ^  =  ttR^  P_/2nR  (8.  1) 

^rad  T 

The  results  are  shown  in  Fig.  8.2a  for  3  atm  and  Fig.  8.2b  for  10  atm. 

(The  5  MW  curve  on  Fig.  8.2b  is  1/10  of  q^ad^*  exactly  the 

radial  heat  flxix  at  the  station,  because  the  radiation  is  emitted  in  all 
directions,  not  only  radially.  However,  the  includsion  of  the  geometrical 
effects  necessitates  a  complicated  calculation  which  is  not  warranted  at 
the  present  stage  of  the  theory.  For  long,  narrow  streamtubes  the  approxi¬ 
mation  is  a  good  one  since  most  of  the  solid  angle  is  radial,  and  the  axial 
energy  transport  out  of  one  station  is  compensated  by  axial  transport  in 
from  other  stations.  For  short,  stubby  streamtubes  like  the  5  MW,  10  atm 
one,  the  approximation  is  not  so  good,  but  these  are  probably  not  desirable 
geometries  to  use. 

The  total  energy  lost  to  the  streamtube  by  plasma  emission  up  to 
the  throat  can  be  obtained  by  integrating  the  curves  of  Fig.  8.  1  from  x  =  0 
to  the  throat  station.  The  average  energy  flux  at  the  streamtube  surface 
is  likewise  obtained  by  integrating  the  curves  of  Fig.  8.2  to  the  throat, 
and  dividing  by  the  throat  station 

^th  ""th 

Q  ,  =  f  TT  R^  P  dx  ,  q  =  — —  J  q  j  dx  (8.  2) 

^rad  J  T  ^rad  x^.  •'  ^rad 

0  0 
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Io=  3.67x10'"  W/m*^ 

Tj  =  12,000  K 

*3.316  kg/m^s 
Pq  *  3.03975  Pa  ( 3 atm) 


0  10  kW 

AlOOkW 
X  5  MW 


/  10  \ 


Throat 


I  X 
X 

Throat 


Throat 


Fig.  8.2a  Station  Radiant  Energy  Flux  of  Laser -Heated 
Streamtube  for  3  atm  Initial  Pressure 


Table  8.  1  presents  these  quantities  for  the  six  cases  calculated, 
as  well  as  the  maximum  value  of  the  energy  flux 


^  (8.3) 

and  the  fraction  of  the  laser  power  radiated.  The  radiative  losses  are 
small  fractions  of  the  power  input,  indicating  that  the  incident  power  is 
converted  efficiently  into  hot  plasma,  which  is  then  available  for  propul¬ 
sion.  The  larger  engines  are  the  most  efficient;  even  though  the  absolute 
amount  of  power  they  radiate  is  larger,  the  fraction  of  the  incident  power 
is  smaller.  The  reverse  is  true  for  the  energy  flux;  the  larger  engines 
have  the  highest  fluxes.  This  is  a  surface  to  volume  ratio  effect  at  any  one 
pressure.  The  higher  fluxes  at  the  higher  pressure  are  a  result  of  the 
shorter  length  of  the  streamtube.  The  magnitude  of  the  fluxes  at  the 
higher  powers  indicate  that  some  form  of  active  heat  protection  will  be 
needed,  either  by  absorbing  some  of  this  energy  in  a  buffer  gas,  or  by 
some  suitable  scheme  of  cooling  the  walls. 

8.3  CONVECTIVE  HEAT  LOSS 

The  hydrogen  boundary  layers  grow  on  the  sides  of  the  core 
streamtube  starting  at  the  initial  station,  with  the  flow  external  to  the 
boundary  layer  given  by  the  core  properties.  For  the  six  cases  calculated 
by  the  core  streamtube  program,  the  Reynolds  number  based  on  conditions 
in  the  core  at  the  throat  station  varies  from  2800  for  the  10  atm,  5  MW  case 
to  8500  for  the  3  atm,  10  kW  case.  This  indicates  that  laminar  flow  will 
prevail,  so  a  laminar  compressible  boundary  layer  model  is  appropriate. 

Since  the  primary  purpose  of  this  calculation  is  to  estimate  the 
heating,  a.nd  to  see  if  a  coupled  calculation  is  necessary,  it  is  appropriate 
to  use  as  simple  a  model  as  will  include  the  important  physical  effects. 

We  have  therefore  assumed  that  a  "local  similarity"  solution  will  be 
satisfactory.  That  is,  at  any  station,  the  boundary  layer  is  calculated  as 
if  it  was  a  similar  boundary  layer  with  the  local  external  properties  as 
parameters.  This  reduces  the  numerical  problem  to  one  of  ordinary, 
rather  than  two  -  variable  partial,  differential  equations,  which  is  a  great 
simplification.  This  approximation  has  proved  to  be  very  useful  in 
calculating  heating  from  high  temperature  gases  ,  and  can  be  found 
discussed  in  books  on  boundary  layer  theory. 

Although  many  cases  of  similar  boundary  layers  have  been  calculated, 
none  known  to  us  is  useful  for  the  present  case  of  very  high  temperature 
hydrogen,  because  they  do  not  include  the  transport  properties  of  hydrogen 
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TABLE  8.  1 


RADIATION  HEAT  LOSSES 


Po 

=  3  atm 

Po 

=  10  atm 

Power  (kW) 

10 

100 

5000 

10 

100 

5000 

Throat  Station  x  , 
th 

(cm) 

3.99 

3.46 

3.24 

0.  54 

0.  51 

0.49 

Power  Loss  Q 

rad 

(kW) 

2.  00 

10.  8 

225 

1.01 

5.  84 

174 

Max  Energy  Flvix 

1.9 

3.6 

11.0 

8.  8 

17.  5 

75 

(qrad^  (kW/cm^) 

max 

Av.  Energy  Flux 

q  -  (kW/cm^) 
rad 

1.  19 

2.  31 

7.  16 

5.  11 

9.95 

43.  1 

%  Power  Loss 

20 

11 

4.  5 

10 

5.  8 

3.  5 

-92- 


in  the  dissociated  and  ionized  state.  We  have  therefore  made  calculations 
including  these  properties,  which  play  an  important  role  in  determining 
the  heating,  since  the  temperature  variation,  and  therefore  the  property 
variation,  across  the  boundary  layer  is  so  large. 


A  simple  set  of  boundary  layer  equations  is  obtained  in  the  standard 
way  by  transforming  from  the  wall  -  oriented  coordinates  s,  y  along  and 
normal  to  the  waU  to  the  similarity  coordinates  I  ,  This  transforma¬ 
tion  is  given  in  terms  of  velocity  u^,  viscosityli  ,  and  density  at  the  edge 
of  the  boundary  layer  by 


I  is)  =  J  P  H  u  ds,  ri(^y)  =  — — - 


(8.4) 


The  dependent  variables  velocity  u  and  total  enthalpy  H  are  written  in 
normalized  form  as 


u/u^  =  f^  ,  H/H^  =  g  (8.  5) 

where  the  subscript  Tl  denotes  differentiation  with  respect  to  'll.  The 
boundary  layer  momentum  and  total  enthalpy  equations  are  transformed 
using  (8.4)  and  (8.  5).  The  "local  similarity"  assumption  is  that  f  and  g 
depend  only  on  'll ;  the  equations  are  then 


)  + 

Tl 


(8.  6a) 


+ 


a 

a-n 


'c  (i-^) 

f  f 

\  Pr  / 

ri  Tiri^ 

(8.  6b) 


where  we  have  defined 

Pr  =  c  li/k  ,  C  =  p  li/p  u  ,  6  =  2d  In  u  /d  In  §  (8.7) 

p  e  e  e 

and  k  denotes  the  thermal  conductivity,  li  the  viscosity. 

The  hypothesis  of  "local  similarity"  has  enabled  us  to  convert  the 
partial  differential  equations  of  momentum  conservation  parallel  to  the  wall, 
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and  energy  conservation,  into  ordinary  differential  equations  in  the  variable 
“n.  The  variable  I ,  which  depends  only  on  distance  along  the  wall,  s,  by 
(8.4),  appears  as  a  parameter  in  u  and  in  8*,  which  involves  the  derivative 
of  Ug  with  respect  to  I  .  The  gas  property  functions  Pr,  C  and  Pg/P  depend 
only  on  the  state  of  the  gas,  so  are  functions  of  p  and  T  (or  p  and  enthalpy 
h  =  H  -  Ug^  / 2)  and  so  can  be  written  as  functions  of  g  at  each  value  of 

^(s).  Thus,  when  the  properties  external  to  the  boundary  layer  are  speci¬ 
fied  as  functions. of  s,  Eq.  (8.  6)  can  be  integrated  to  find  f('n),  g('n). 

2 

Two  further  simplifications  are  useful.  Since  Ug  /Hg  is  small,  and 
Pr  is  near  unity  for  hydrogen  over  the  range  of  interest,  we  will  put  Pr  =  1 
in  (8.  6b).  Further,  the  third  term  in  (8.  6a),  the  pressure  gradient  term, 
is  small  so  we  will  approximate  Pg/P  by  h/h  .  When  this  is  expressed  in 
terms  of  g,  a  new  pressure  gradient  parameter  P  appears. 


d  Inu 


_ e_  ■ 

d  In  5 


Now  the  boundary  layer  equations  are 


+  +P  (g  -  f^  )  =  0 


T1 


r\r\ 


r] 


(8.8) 


(8.9a) 


(Cg  )  +  fg  -  0  (8.9b) 

T) 

The  boundary  conditions  are  no  slip  at  the  wall,  a  given  wall 
enthalpy  h^  based  on  the  wall  temperature  T^,  and  approach  to  the  exter¬ 
nal  conditions  at  large  T) ; 


(0)  =  0,  g  (0)  =  h^H^.=  g^.  f^  (»)  =  g  (CO)  =  1  (8.  10a) 

For  the  case  where  there  is  no  flow  through  the  wall,  we  have  the  fifth 
condition  that  the  normal  velocity  is  zero  there,  which  is  given  by 

f  (0)  =  0  (8.  10b) 
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When  the  external  conditions  pg,  (jig,  Ug  are  given  from  the  core 
streamtube  solutions,  and  the  p  and  \x  are  defined  in  terms  of  h,  then 
(8.  9)  and  (8.  10)  can  be  integrated.  The  p  (h)  function  was  found  from  the 
equilibrium  properties  for  hydrogen  given  in  Ref.  33,  as  was  the  function 
h(T).  The  viscosity  was  found  from  Ref.  32  as  a  function  of  T  and  p, 
completing  the  necessary  information.  The  integration  was  performed  by 
the  method  of  quasi -linearization. 


The  results  of  the  integration  for  Tg  =  14,  OOOK  p  =  3  atm,  are 
shown  in  Fig.  8.  3  as  plots  of  the  shear  and  heat  transfer  parameters 
f.^.j.|{0)  and  g.q  (0)  against  the  pressure  gradient  parameter  p,  for  various 
values  of  g.^.  These  are  related  to  the  shear  stress  and  heat  flux  q^ 
by 


T  = 
W 


2 

p  |i  u 
w  w  e 


JU 


w 


p  li  u  H 
w  w  e  e 

<j2lPr 


gr^  (0) 


(8.11) 


In  addition  to  the  solid  curves  for  equilibrium  hydrogen  properties, 
the  dashed  curves  show  the  results  of  assuming  that  the  p  jj  product  does 
not  vary  across  the  layer,  ie.,  C  =  l.  For  the  small  values  of  wall  to 
external  enthalpy  ratio  gw  in  which  we  are  interested  (g^  ^  0.  02  for 
Tw  =  1, 300  K,  the  melting  point  of  copper)  it  is  very  important  to  consider 
the  variation  of  C.  This  is  not  surprising,  since  for  3  atm,  Tg  =  14,  OOOK, 
Tw  =  1. 300K,  C  varies  from  unity  (by  definition)  at  the  outer  edge  of  the 
boundary  layer,  to  7.4  at  the  wall.  Another  notable  feature  is  the  very 
slight  dependence  of  g.^  (0)  on  pressure  gradient  parameter  p,  especially 
for  g.^  ~  0.  02 . 


Based  on  calculations  like  those  shown  in  Fig.  8.  3,  and  the  use  of 
Eq,  (8.  11),  the  heat  transfer  distribution  along  a  1, 300 K  wall  in  contact 
with  the  plasma  core  was  calculated  for  the  3  atm,  10  kW  case.  Because 
conditions  in  the  core  for  the  100  kW  and  5000  kW  cases  are  only  slightly 
different,  the  resulting  heating  rate  is  also  a  good  estimate  for  these  two 
cases  at  3  atm.  The  distribution  of  q.^  is  shown  in  Fig.  8.4,  and  gives 
fluxes  from  10  kW/cm^  to  a  peak  of  34  kW/cm^  near  the  throat.  If  these 
values  are  integrated  over  the  surface,  the  total  convective  heat  loss  rate 
is  found  from 


2Tr  R  q  dx 
w 


(8.12) 


The  values  that  result  for  the  10  kW,  100  kW  and  5000  kW  cases  (each  of 
which  have  a  different  radial  size),  for  integration  from  x  =  1. 5  to  x^^* 
17,53  and  374  kW  respectively. 
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Fig.  8,3a  Shear  Parameter  for  Hydrogen  Boundary  Layer 
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Fig,  8.3b  Heat  Transfer  Parameter  for  Hydrogen  Boundary  Layer 
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It  is  iirtmediately  evident  that  this  uncoupled  calculation  produce 
unrealistic  results  for  both  the  10  kW  and  100  kW  cases,  since  the  calculated 
heat  loss  rate  is  a  large  fraction  of  the  total  power.  Such  a  loss  rate  would 
decrease  the  core  temperature,  which  in  turn  would  decrease  the  loss. 

For  the  low  power  cases,  then,  the  convective  loss  must  be  coupled  to  the 
core  calculation  for  a  realistic  model  of  the  fluid  mechanics.  For  the 
500  kW  case,  however,  the  loss  calculated  is  only  7.  5%  of  the  laser  power, 
which  is  probably  not  a  large  enough  loss  to  require  coupling  to  the  core. 

F or  this  case,  the  heating  distribution  shown  in  Fig.  8.  4  is  probably  realistic. 
This  certainly  indicates  the  need  for  some  protection  of  the  rocket  walls, 
since  rates  of  30  kW/cm^  cannot  be  handled  by  simple  means.  However, 
the  larger  rates  do  occur  near  the  throat,  where  the  streamtube  is  of  small 
radius,  so  the  area  to  be  protected  is  not  large.  The  throat  radius  for  the 
3  atm,  5000  kW  case  is  about  1  cm. 

Realistic  heat  fluxes  for  the  lower  power  engines  can  only  be 
determined  by  a  coupled  calculation  of  the  core  and  boundary  layer,  which 
was  not  undertaken  under  the  present  contract. 

The  boundary  layers  for  the  10  atm  cases  were  not  calculated,  but 
for  laminar  heating,  the  heat  fluxes  should  scale  like  the  square  root  of 
pressure,  increasing  them  by  a  factor  of  about  2.  But  the  areas  of  the 
streamtube  surface  are  nearly  10  times  smaller  because  of  the  short 
length,  so  the  values  of  Q^.  should  be  about  one  fifth  of  those  for  3  atm. 

This  means  the  coupled  calculation  is  still  necessary  for  the  10  kW  case 
at  10  atm,  might  not  be  necessary  for  the  100  kW  case,  and  of  course  is 
not  necessary  for  the  5000  kW  case.  The  heat  protection  of  the  walls 
will  be  even  more  important  though,  since  q^  is  higher  by  a  factor  of  2. 

One  method  of  heat  protection  is  the  injection  of  cold  gas  normal  to 
the  wall,  which  alters  the  temperature  profile  and  blocks  some  of  the 
heat.  Some  sample  calculations  have  been  made  to  illustrate  the  effect 
of  this  method  of  reducing  the  heat  flux.  The  only  change  in  the  calcula¬ 
tion  is  the  introduction  of  a  non-zero  value  of  f  (0)  to  replace  the  boundary 
condition  (8.  10b),  because  now  the  mass  flux  through  the  wall  is  not  zero, 
but  is  related  to  f  (0)  by 


(pv)w  =  (0)/V2? 


(8. 13) 


This  introduces  another  parameter  into  the  equations,  f  (0).  The  results 
of  two  calculations,  at  stations  x  =  1  and  x  =  2.  7,  are  given  in  Fig.  8.  5  as 
curves  of  g  (0)  against  f  (0).  The  left  end  of  these  curves  is  the  zero 
injection  or^olid  wall  case  f  (0)  =  0.  It  can  be  seen  that  the  heat  transfer 
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Effect  of  Normal  Hydrogen  Injection  on  Heat 
Transfer  Parameter  for  Hydrogen  Boundary  Layer 


parameter  (0)  drops  sharply  as  f  (0)  decreases,  which  means  the 
injection  velocity  increases  by  (8.  13).  If  we  look  at  injection  with 
f  (0)  =  -  1,  (0)  at  X  =  1  is  reduced  by  a  factor  of  4,  leading  to  a  heat 

flux  of  about  2,5  kW/cm^  for  an  injection  rate  of  C,014  gm/cm^-s, 
from  (8,  13),  At  x  =  2.  7,  g^  (0)  is  reduced  by  a  factor  of  2,  yielding 
17  kW/cm^  with  an  injection  rate  of  0.017  g/cm^-s. 

To  see  the  magnitude  of  this  rate  of  injection,  one  may  compare 
it  to  the  mass  flow  in  the  core,  which  is  0.  3316  g/cm^-s  over  the  initial 
area.  Thus  for  the  10  kW  case,  the  core  mass  flow  is  9.  OE-3  g/s,  while 
an  injection  of  0.  015  g/cm^-s  over  the  surface  area  of  the  streamtube 
up  to  the  throat  is  22E-3  g/s,  more  than  twice  the  core  mass  flow.  On 
the  other  hand,  for  the  5000  kW  case,  the  core  mass  flow  is  500  times 
larger,  4.  5  g/s,  while  the  area  is  only  ^  500  =  22  times  larger,  so  the 
injected  mass  would  be  about  0,49  g/s,  or  only  11%  of  the  core  mass 
flow.  The  surface  to  volume  effect  favors  injection  in  the  large  engines, 
and  it  appears  to  be  one  possibility  for  reducing  the  wall  heat  load  by  a 
factor  of  2  or  more.  Notice  also,  that  large  reductions  in  heat  flux  are 
only  needed  in  the  limited  regions  of  high  heating,  not  over  the  whole 
streamtube  surface,  so  injection  is  not  needed  over  the  whole  nozzle 
surface,  which  reduces  the  mass  which  must  be  injected. 

The  purpose  of  the  present  calculation  was  to  show  how  much 
heating  reduction  can  be  achieved  by  mass  injection.  Considerable  addi¬ 
tional  calculation  would' be  needed  for  an  engineering  analysis  of  the 
proper  amount  and  distribution  of  injection  for  efficient  heat  protection. 

8.4  ABSORPTION  BY  PARTICLES  IN  BUFFER  GAS 

It  appears  that  for  specific  device  configurations  the  radiation 
escaping  from  the  plasma  will  be  of  sufficiently  high  intensity  that  some 
sort  of  shielding  mechanism  will  be  required  in  order  to  maintain  wall 
integrity.  One  possibility  would  be  to  introduce  small  carbon  particles 
into  a  buffer  gas  flowing  around  the  plasma.  In  this  event  the  carbon 
particles  would  intercept  the  plasma  radiation,  heat  up  and  conduct  heat 
energy  to  the  surrounding  buffer  gas.  These  particles  would  also  radiate 
with  an  intensity  corresponding  to  their  own  characteristic  temperature 
and  could  further  reach  vaporization  temperatures  depending  on  the 
particle  flow  time  and  their  radiative  flux  loading. 

In  the  following  analysis  first  order  estimates  are  provided  for 
the  particle  mass  loading  required  to  protect  the  walls  against  plasma 
radiation.  The  particles  are  assumed  to  have  constant  absorption  proper - 
ties,  velocities  and  size.  An  understanding  of  the  phenomenology  of  the 
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plasma /particle /buffer  gas  interaction,  including  the  all-important 
question  of  particle  survival,  will  require  a  significantly  more 
sophisticated  analysis  of  the  mixing  region  between  the  plasma  and 
buffer  gas  flows.  When  such  an  analysis  is  performed  the  required 
particle  mass  loadings  can  be  more  accurately  defined. 

Small  carbon  particles,  of  the  order  of  1  ijm,  are  very  efficient 
absorbers  of  radiation  in  the  visible  and  near  infrared  wavelength 
regions,  approximately  acting  as  black -body  absorbers.  A  large  fraction 
of  the  hydrogen  plasma  radiation  will  unfortunately  originate  in  the  far  UV 
where  little  is  known  about  the  absorptive  properties  of  carbon.  In  the 
following  analysis  it  will  be  assumed  that  the  carbon  particles  absorb  as 
black  bodies  at  all  wavelengths. 

The  configuration  considered  is  a  cylindrical  plasma  of  radius  ri 
surrounded  by  a  particle  seeded  buffer  flow  of  thickness  /\r.  The 
particles  are  taken  to  be  spherical  with  radius  Rp  and  to  have  a  constant 
number  density  np  in  the  buffer  flow.  Each  particle  has  an  absorption 
cross  section  of  ttR  ^  and  thus  the  optical  depth,  t,  across  the  buffer 
layer  is  defined  by 


T  =  n  ttR  ^ 

P  P 

An  optical  depth  of  3  is  sufficient  to  absorb  95%  of  the  plasma  radiation 
and  this  condition  is  met  when  the  particle  number  density  is  given  by 


n  =  - r - 

P  TTRp  Ar 


R  ^  At 
P 


(8.  15) 


It  is  perhaps  more  interesting  to  examine  this  quantity  in  terms 

of  the  total  mass/cc  of  particles  in  the  buffer  flow,  M  , 

P 


M  =n  p  (4/3)ttR„' 
P  P  P  P 


(8.  16) 


where  p  is  the  specific  density  of  carbon  which  is  ~2  gm/cc.  Combining 
Eq.  (8.  1%)  and  (8.  16)  results  in 


R  R 

M  =  4  pp  ~  8  — ^ 

P  Ar  Ar 


(8.17) 
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Equation  (8.  17)  can  readily  be  extended  to  deduce  the  mass  flow 
of  particles  required  to  absorb  the  plasma  radiation. 


where 


M  =  M  u^  A, 
P  P  h  b 


=  TT  (2  Ar  +  (Ar)  ) 


(8.  18) 

(8.  19) 


and  ub  is  the  buffer  flow  velocity.  Since  it  is  desirable  to  maintain  a  con¬ 
figuration  such  that  Ar/i’i  «  1>  the  approximation 


2  Trr^  Ar 

is  made.  Combining  Eq.  (8.17),  (8.18)  and  (8.20)  results  in 


(8.20) 


It  can  be  seen  from  Eq.  (8.  21)  that  as  long  as  A  ^  «  ^1>  the 
particle  mass  flow  required  to  absorb  the  plasma  radiation  is  independent 
of  the  buffer  flow  thickness  and  linearly  proportional  to  both  the  particle 
and  plasma  radii.  In  the  interest  of  optimizing  thrust  one  wishes  to 
minimize  the  buffer  flow;  hence,  it  is  clear  that  the  smallest  particle 
size  consistent  with  the  assumption  of  blackbody  absorption  should  be  used. 

Another  interesting  feature  is  that  since  the  required  buffer  flow 
is  linearly  proportional  to  plasma  radius,  the  ratio  of  buffer  flow  to  plasma 
flow  is  size  dependent;  i.  e. 


8  p  R  u, 
V^p  p  b 

P  1  ^1  1 

^pl  1  pi 


(8.  22) 


where  u  is  the  plasma  flow  velocity.  For  example,  for  a  14,  000°K  hydrogen 
plasma  eR:  a  pressure  of  3  atm. 


10 


„  R  u, 
7  p  b 


r.  u  , 
1  pi 


(8.23) 
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Note  that  relationship  (8.23)  implies  that  for  small  devices,  i.  e. 
r  j  <  1  cm,  the  buffer  mass  loading  would  of  necessity  be  orders  of 
magnitude  larger  than  that  of  the  plasma  for  reasonable  particle  sizes 
unless  ub«u  Thus  these  preliminary  numbers  suggest  that  particulate 
absorption  may  not  represent  a  viable  technique  for  wall  protection  in  the 
flowing  plasma  configuration.  Note  that  even  if  the  buffer  flow  had  a 
velocity  two  orders  of  magnitude  lower  than  the  plasma,  a  condition 
difficult  to  maintain  along  the  flow  length,  there  may  still  be  an  unacceptable 
constraint  on  particle  size. 

\ 

Although  the  analysis  presented  above  is  admittedly  crude  it  would 
suggest  that  particle  shielding  would  provide  too  large  a  weight  penalty  for 
proposed  operating  conditions.  A  more  detailed  analysis  would  be  required 
to  confirm  this.  Furthermore  the  possibility  of  operating  under  different 
conditions,  for  instance  at  higher  plasma  pressures,  should  also  be 
considered. 
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9.  PERFORMANCE  OF  FLOWING  CORE  ROCKET 


The  performance  of  the  flowing-core  laser -heated  rocket  is  affected 
by  many  more  variables  than  the  ones  considered  in  this  report.  Most  of  the 
effort  reported  here  has  dealt  with  the  fluid  mechanics,  laser  radiation 
absorption  and  plasma  emission  of  the  flowing  core,  and  the  convective  heat 
transfer  from  the  core  streamtube.  Based  on  this  work,  only  the  performance 
obtained  from  the  core  streamtube  can  be  obtained,  that  is,  the  thrust  and 
specific  impulse  of  the  core  gas  when  it  is  expanded  to  some  exit  condition. 

The  resulting  values  will  be  high  estimates  because  no  losses  past  the  throat 
are  considered,  nor  has  the  convective  heating  loss  been  properly  coupled  to 
the  core  streamtube  for  the  lower  power  levels  as  yet.  On  the  other  hand, 
added  thrust  which  might  be  supplied  by  a  buffer  gas  layer  used  for  heat 
protection  is  not  included  either.  With  the  understanding  of  the  partial  nature 
of  the  performance  which  can  be  calculated  here,  we  can  describe  the  method 
for  the  calculations. 


We  begin  at  the  throat  of  the  streamtube  with  the  enthalpy  h^j^  and 
velocity  u^j^  found  from  the  program  described  in  the  section  on  the  Flowing 
Core  Streamtube.  The  flow  is  then  expanded  isentropically  to  a  specified 
exhaust  pressure  Pg^  by  linear  interpolation  in  hydrogen  property  tables  in 
Ref.  33.  This  provides  the  exhaust  conditions  of  enthalpy  and  temperature 
hex*  Tgjj.  Conservation  of  total  enthalpy  during  the  expansion  (no  losses) 
then  provides  the  exhaust  velocity  as 


u 

ex 


-  h  ) 
ex 


The  thrust  is  then 


Th 


(p  -  p  )  A  +  m  u 
ex  am  ex 


ex 


(9.1) 


where  p  is  the  ambient  pressure.  We  will  give  the  so-called  ideal  thrust, 
■^axn  ^ 

for  which  Pg^  =  Pam»  that 


Th  =  mu  ,  I  =  u 

ex  sp  ex 

give  the  thrust  and  specific  impulse,  where  g 


/g  (9.2) 

is  the  acceleration  of  gravity. 
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The  conditions  at  the  throat,  as  calculated  by  the  core  streamtube 
program,  the  mass  flow  rate  from  Table  7.2  and  the  entropy  s  from 
Ref.  33  are  given  in  Table  9.  1.  The  exhaust  conditions,  Igp  and  Th  for 
Pex  =  1  atm,  as  found  from  Ref.  33,  and  Eqs.  (9.  1)  and  (9.^)  are  given 
in  Table  9.  2.  The  maximum  thrust  is  obtained  by  expanding  to  pg^  =  0, 
since  all  energy  is  then  converted  into  velocity,  hg^j.  being  zero.  This 
limiting  case  is  given  in  Table  9.3,  as  obtained  from  Eqs.  (9.  1)  and  (9.2) 
by  setting  hgx  =  0. 

Table  9.2  shows  that  the  expansion  to  1  atm  of  the  3  atm  cases  is 
really  a  very  small  expansion,  since  the  throat  is  at  about  1.  7  atm.  There¬ 
fore,  the  gas  in  the  exhaust  is  still  very  hot  (near  13,000  K)  and  its 
enthalpy  is  only  a  few  percent  less  than  at  the  throat.  It  is  this  few  per¬ 
cent  which  is  converted  to  additional  velocity,  so  the  increase  from  the 
throat  velocity  to  the  exhaust  velocity  is  also  only  a  few  percent.  However, 
this  velocity  is  sufficient  to  produce  Igp  in  the  range  of  1400  to  1500  s. 
Because  of  the  small  mass  flow  rate,  the  thrust  is  in  the  range  0.  1  to  66  N 
(0.03  to  15  lbs.). 

The  10  atm  cases  have  a  greater  expansion  ratio,  since  the  throat 
is  at  5.6  to  5.  8  atm,  .so  they  have  a  greater  conversion  of  enthalpy  to 
velocity,  exhausting  somewhat  cooler  gas,  at  speeds  nearly  double  the 
throat  speed.  They  then  produce  higher  values  of  Ig  ,  in  the  range  2400 
to  2500  s,  and  somewhat  larger  thrusts,  0.2  to  100  N  (0.05  to  25  lbs.  ). 

For  the  maximum  possible  expansion,  in  Table  9.  3,  we  see  much 
less  difference  between  the  3  atm  and  10  atm  cases,  since  the  expansion 
ratio  is  infinite  for  both.  They  all  yield  Igp  in  the  range  3000  to  3300  s. 

In  fact,  the  100  and  5000  kW  cases  have  nearly  the  same  at  Igp  at  both 
pressures.  Only  the  highest  percent  loss  case,  10  kW,  has  noticeable 
difference  in  Igp  with  pressure.  The  thrusts  for  the  vacuum  exhaust 
cases  are  also  higher  than  the  atmospheric  exhaust  cases  since  thrust 
is  proportional  to 

Even  though  a  number  of  loss  mechanisms  were  omitted  from 
this  calculation,  and  they  will  reduce  the  performance  calculated  here,  it 
seems  quite  clear  that  values  of  Igp  above  1000  s  can  be  obtained,  even 
over  a  3  atm  to  1  atm  expansion. 

The  thrusts  produced  may  seem  rather  small,  because  of  the  very 
small  mass  flow  rates.  However,  the  Th/P  ratio  for  the  vacuum  expansion 
core  is  in  the  range  20  to  30  mN/kW,  which  seems  to  be  typical  for  electric 
propulsion  schemes,  to  which  the  laser-heated  rocket  is  the  most  comparable. 
Of  course,  the  laser -heated  rocket  does  not  have  to  carry  its  own  power 
supply,  so  its  weight  for  a  given  mission  should  be  smaller  than  an  electric 
propulsion  engine. 
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THROAT  CONDITIONS 
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PERFORMANCE  FOR  =  0  (VACUUM) 
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10.  STATIONARY  CORE  ENGINE 


10.  1  INTRODUCTION 

An  alternative  to  the  flowing  plasma  core  type  of  rocket  is  one 
borrowed  from  the  concept  of  a  gas  -core  nuclear  rocket.  In  that  type 
of  device  a  gaseous  core  of  fissionable  material  is  maintained  in  a 
stationary  region  in  the  center  of  the  rocket  while  fissioning  at  high 
temperatures.  It  is  surrounded  by  a  flowing  propellant  gas  which  confines 
the  reacting  core,  and  is  heated  by  thermal  radiation  and  convective 
processes.  The  nuclear  fuel  burns  above  100,  OOOK,  and  it  must  be  kept 
at  least  at  atmospheric  density  to  maintain  criticality,  so  it  operates  at  a 
core  pressure  of  several  hundred  atmospheres. 

An  analogous  laser-heated  rocket  would  have  a  stationary  plasma 
core,  which  absorbs  the  laser  radiation  and  is  thereby  maintained  at  high 
temperature.  In  turn,  this  core  re -radiates  the  absorbed  laser  energy 
to  the  surrounding  propellant  gas,  which  serves  to  confine  the  core  and 
produce  thrust,  just  as  in  the  gas -core  nuclear  rocket.  This  concept  may 
be  compared  with  both  the  gas -core  nuclear  rocket  and  the  flowing  core  laser- 
heated  rocket. 

Compared  to  the  former,  there  is  no  criticality  requirement 
imposed  on  the  core.  However,  laser  energy  absorption,  which  is  the 
mechanism  for  putting  energy  into  the  rocket,  is  dependent  on  the  core 
density.  The  more  dense  the  core,  the  better  it  absorbs,  and  the  shorter 
it  can  be  in  the  direction  of  the  laser  beam.  But  there  is  more  freedom  in 
the  range  of  densities,  and  therefore  pressures,  in  the  core  for  the  laser - 
heated  case  than  for  the  gas -core  case.  Another  difference  is  in  the 
amount  of  core  plasma  which  can  be  allowed  to  flow  out  the  rocket  nozzle. 

For  nuclear  rockets,  one  would  like  to  have  as  little  core  material  as 
possible  entrained  by  the  propellant  gas  and  carried  out  the  nozzle,  since 
it  represents  a  loss  of  energy  of  fission  as  well  as  of  radio-active  material. 

In  the  laser  rocket,  ’’detrainment”  of  core  material  is  not  so  serious, 
although  it  does  represent  a  loss  of  propulsive  efficiency,  since  it  will 
leave  the  rocket  at  rather  high  temperatures  thus  wasting  energy.  However, 
if  we  let  too  large  a  fraction  of  the  core  flow  out  with  the  propellant  gas,  we 
are  back  to  the  flowing  core  concept.  Preventing  this  detainment  is  not  so 
easy  in  the  laser  rocket,  since  the  core  may  not  be  of  high  density  relative 
to  the  propellant  gas,  as  in  the  nuclear  core. 
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Compared  to  the  flowing -core  laser -heated  rocket,  there  are 
several  new  physical  mechanisms  which  enter.  One  is  the  detrainment 
problem  referred  to  above,  which  acts  against  the  confinement  of  the  core. 
The  densities  of  the  core  and  propellant  may  be  nearly  equal,  which  makes 
confinement  difficult.  Further,  the  heating  of  the  propellant  and  its  sub¬ 
sequent  density  decrease  occur  while  it  is  flowing  past  the  core,  making 
detrainment  easier  as  it  flows.  A  second  important  mechanism  is  the 
efficient  emission  of  radiation  by  the  core,  and  its  efficient  absorption  by 
the  cooler  propellant  gas.  In  the  flowong-core  scheme,  emission  by  the 
core  is  to  be  minimized  as  a  loss  of  energy  from  the  primary  propellant, 
while  in  the  stationary -core  this  emission  is  the  main  process  by  which 
energy  is  input  to  the  propellant.  Likewise,  absorption  by  the  propellant 
gas  is  now  a  primary  source  of  propulsion  rather  than  a  necessary  heat 
protection  device,  as  in  the  flowing -core  rocket.  In  fact,  one  of  the  main 
complications  of  the  stationary -core  rocket  is  in  the  control  of  the  flowing 
propellant  gas.  This  is  achieved  primarily  by  control  of  the  flow  area 
transverse  to  its  direction  of  motion,  i.  e.,  the  cross-section  area  of 
the  propellant  gas  streamtube.  The  same  is  true,  of  course,  in  the  flowing- 
core  case;  its  flow  properties  are  largely  controlled  by  the  cross-sectional 
area.  But  there  is  one  crucial  difference.  In  the  flowing -core  case,  the 
propellant  gas  absorbs  energy  from  a  laser  beam  propagating  in  the  axial 
direction,  so  its  absorption  properties  depend  on  its  axial  size.  Thus  the 
control  of  flow  and  absorption  depend  largely  on  two  independent  geometrical 
parameters,  the  cross-section  and  the  axial  length.  However,  in  the 
stationary-core  rocket,  the  propellant  gas  absorbs  its  radiant  energy 
mainly  from  emission  in  the  transverse  direction,  the  same  direction  which 
controls  the  flow.  This  means  that  flow  control  will  be  intimately  connected 
with  absorption,  and  the  propellant  gas  flow  area  must  be  designed  with 
regard  for  both  the  aixial  pressure,  velocity  and  density  distributions  and 
for  the  gas  depth  for  efficient  absorption. 

It  appears,  then,  that  the  analysis  of  the  stationary -core  rocket  is 
considerably  more  complex  than  that  of  the  flowing-core  rocket,  with 
increased  emphasis  on  plasma  emission  and  subsequent  absorption,  and 
a  more  complicated  interaction  between  flow  and  radiation. 

10.2  ORDER  OF  MAGNITUDE  ANALYSIS  ' 

The  quantitative  evaluation  of  a  stationary-core  rocket  can  be 
started  by  using  some  simple  relations  among  the  geometric  and  physical 
parameters  which  serve  to  indicate  the  physics  which  goes  into  such  an 
analysis,  as  well  as  obtain  some  order -of-magnitude  estimates  of  the 
parameters . 
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Taking  a  cylindrical  core  of  radius  R^,  and  length  L,  the  laser 
power  P  and  intensity  I  are  related  by 


P  =  TT  R  ^  I  (10.  1) 

c 

The  absorption  of  laser  energy  occurs  with  absorption  coefficient  a 

function  of  temperature  T.  The  volume  emission  by  the  plasma, 
also  depends  on  T.  In  steady  state  the  plasma  must  emit  the  full  power 
P,  so  that 

P  =  TTR  ^  LP„  (10.2) 

c  ^ 


Equations  (10.  1)  and  (10,2)  give 

L  =  I/P  (10.  3) 

c 

The  absorption  and  emission  in  the  core  are  related  by  the  steady  state 
requirement  that  the  absorption  kj^I  equals  the  emission  P^: 


k,  I  =  P  (10.  4) 

Ij  c 

Equations  (10.  4)  and  (10.  3)  show  that  with  these  estimates  ,  the  core 
length  L  is  equal  to  the  absorption  coefficient  k^. 

Of  the  emitted  radiation  per  unit  volume  P^,  some  portion  P^<Pc 
will  be  absorbed  by  the  propellant  gas,  depending  on  the  depth  of  this  gas, 
and  its  absorption  coefficient  kp.  If  mp  is  the  mass  flow  rate  of  the 
propellant  gas  ,  the  energy  per  unit  mass  added  is  tt  R^,^  •  If  its 

initial  enthalpy  is  hi,  then  its  heated  total  enthalpy  is 

H  h.  +  u?/2  +  ttR^  LP^'Vm^  (10.  5) 

h  1  I  c  c  P 


This  gas  is  then  expanded  at  constant  entropy  to  the  chosen  exhaust  pressure, 
p^^,  yielding  an  exhaust  enthalpy  and  an  exhaust  velocity 


u 

ex 


and  associated  specific  impulse 


(Hh  - 
^ex/g  =  fgp- 


(10.6) 
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With  these  relations,  estimates  of  the  parameters  of  a  rocket 
can  be  made.  Suppose  we  give  the  power  P,  core  temperature  Tc,  core 
pressure  p^,  exhaust  pressure  pex  and  specific  impulse  desired,  which 
defines  Uex.  From  Tc  and  the  radiation  properties  of  hydrogen  we  can 
find  Pj,  and  kj^.  Then  (10.4)  gives  I,  (10.3)  give  L  and  (10.  1)  give  Rc, 
so  the  core  size  is  determined. 

The  mass  flow  rate  of  propellant  mp  may  now  be  found  from  Ug^. 
Equation  (10.6)  gives  -  hex*  If  we  assume  uj^/s  is  a  small  contribution 
to  this  enthalpy  difference,  then  we  know  the  desired  difference  between 
the  enthalpy  of  the  propellant  gas  at  the  end  of  the  heating  zone,  which  we 
will  take  to  be  at  constant  pressure  Pc,  and  its  enthalpy  at  the  exhaust 
pressure  Pg^*  These  two  states  are  connected  by  constant  entropy.  The 
conditions  uniquely  define  both  states,  so  the  values  of  Tj^,  hjj,  Tg^^,  hgjj. 
can  be  found.  Again  assuming  the  velocity  contribution  to  the  total  enthalpy 
is  small,  we  take  hj^  -hj  to  be  energy  added  to  the  gas.  Now  we  estimate 
the  volumetric  absorption  P^  as  a  fraction  of  the  emission  Pq.  Then 
(10.  5)  yields  the  propellant  rnass  flow  as 

.  mp  =  L  P';/(h^-h.  )  (10.7) 


where  the  initial  enthalpy  is  assumed  known.  The  ideal  thrust  can  also  be 
found  now  as 


Th  =  m,^  u  (10.  8) 

P  ex 

To  check  the  neglect  of  the  velocity  compared  to  the  enthalpy  during 
heating,  and  to  estimate  detainment  effects,  we  need  to  define  the  area  in 
which  the  propellant  gas  flows.  This  is  in  turn  related  to  the  absorption 
length  of  plasma  radiation  in  the  gas,  and  so  depends  on  the  detailed  proper¬ 
ties  of  the  gas.  If  the  radius  of  the  outer  edge  of  the  propellant  gas  is  Rp 
at  any  station,  its  velocity  is 


Up  =  liip/pp  tt  (Rp  -Rp^  )  (10.9) 

and  pp  is  known  from  the  given  gas  state  at  the  inlet,  and  the  end  of  heating 
state,  which  was  found  above. 

The  detrainment  occurs  in  the  shear  layer  between  the  stationary  core 
and  the  flowing  propellant,  where  the  two  gases  mix.  If  the  thickness  of  the 
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mixing  layer  on  the  propellant  gas  side  is  6,  then  the  core  gas  mass  flow 
there  is 


m 

c 


2TrR 

c 


P'  u'  dy 


(10.  10) 


I  I 

where  p  ,  u  are  the  density  and  velocity  of  the  core  gas  in  the  mixing 
layer.  By  normalizing,  this  can  be  written 


m  -  2TT  R  S  p  a  f  - (10.11) 

c  c  c  Pj  P  u  \  5  / 

0 

The  factor  in  front  gives  the  order  of  magnitude  of  the  mass  flow  rate  out 
of  the  core  at  any  station,  and  the  total  core  detrainment. 

The  thickness  6  is  determined  by  a  balance  between  the  diffusion 
of  core  gas  into  propellant  gas  and  the  convection  of  the  propellant  past  the 
core.  This  can  be  estimated  by  combining  the  static  diffusion  length //dT 
with  the  time  it  takes  for  the  propellant  gas  to  flow  past  the  core, 

t  =  ^ dx/up.  Here  D  is  the  diffusion  coefficient  of  core  gas  into 
propellant  gas.  Thus  at  the  end  of  the  core 


«  =  dx/Upj^^^  (10.12) 

'o' 

To  evaluate  the  integral,  we  will  use  a  linear  increase  of  up  from  Up.  to 
uph,  so 


a 

V“ph“‘^Pi 


In 


‘^Pi/ 


(10.  13) 


Now  the  factor  in  frcxit  of  the  integral  in  Eq.  (10.  11)  can  be  calculated  as 
an  estimate  of  detrainment. 


10.3  NUMERICAL  EXAMPLE 

As  an  example  of  the  use  of  these  formulas  for  estimating  the 
performance  of  a  small  stationary  core  engine  we  may  take  as  a  core: 


-115- 


P  =  10  kW,  p.  =  P  =  3  atm,  T  =  14,  000  and  18,  OOOK 
1C  c 

p  =1  atm,  =  1000  s 

sp 


Then  the  core  conditions  are: 


T^(K) 

14,  000 

18,000 

p^  (kW/cm^) 

56.4 

200 

(cm"S 

1 

3 

I  (kW/cm^) 

56 

67 

(cm) 

.24 

.22 

L  (cm) 

1 

.  33 

Pg  (g/cm^) 

2.  33  E-6 

1.25  E-6 

p*  (kW/cm^)- 
c 

29 

90 

TTR2g  L  P*  (kW) 

5.2 

4.  5 

The  emitted  radiation  which  is  assumed  to  be  absorbed  by  the 
propellant  gas,  P^,  is  that  below  0.  1  ijm.  The  remaining  energy  is  lost. 
This  is  a  rough  estimate  based  on  our  present  knowledge  of  absorption 
by  hydrogen,  and  ignores  IR  absorption. 

Turning  to  the  propellant  gas,  Igp  =  1000  s  corresponds  to 
Uex  =  980E5  cm/s  and  a  kinetic  energy  per  unit  mass  of  4.  8E7J/kg. 
Hydrogen  properties  at  pgx  =  1  atm  then  show  that  at  Tg^  =  5000  K, 

Ilex  =  0-  93  E8  J/kg  at  an  entropy  of  Sg^  =170  J/kgK,  which  leads  to 
Hex  =  1.41  e*  J/kg.  The  corresponding  point  at  p^  =  3  atm  is  at  7000 K, 
hh  =  1. 43E8  J/kg,  Sh=  170  J/kgK,  pj^=5.  3E-6  g/cm^,  which  gives 
the  state  of  the  propellant  gas  after  it  is  heated.  Knowing  the  power  to 
be  absorbed  to  heat  from  =  BOOK,  =  2.  10  E8  J/kg,  =  3  atm  to  this 
state,  we  find  the  mass  flow  rate  of  propellant  gas,  mp,  as  0.015  and 
0.013  g/s  for  the  two  core  temperatures. 

To  find  the  propellant  gas  speeds  and  core  detra inment  we  must 
now  make  some  assumption  about  absorption  depths  to  absorb  the  UV  radia¬ 
tion  A  reasonable  estimate  is  0.  1  cm,  leading  to  Rp  =  +  0*  1* 
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F or  this  constant  area  flow  of  propellant  gas,  the  value  of  up  p  comes 
from  Eq,  (10.  9),  and  since  the  density  is  known  from  the  constant 
pressure  and  known  temperatures  Ti  =  300K,  =  7000K,  the  velocities 

can  be  found,  as  can  the  constant  factor  in  Eq,  (10.  11),  us ing  (10.  13) 
for  6.  The  diffusion  coefficient  D  is  evaluated  at  =  7000  K  by  scaling 
the  value  for  H,  diffusion  Ln  Reference  32  down  to  3  atm. 


T  (K) 
c 

14,000 

18,000 

mp(g/s) 

1.5E-2 

1.3E-2 

Pp  Up(g/cm^^  s) 

8.2E-2 

7.7E-2 

Up^(cm/  s) 

332 

312 

Uph'^m/s) 

15.5E3 

14.5E3 

D(cm^/  s) 

230 

230 

6  (cm) 

0.23 

0.14 

m  ~  2iTR^6p^Upj^(g/s) 

1.2E-2 

3.5E-3 

Core  Mass  (g) 

2.9E-7 

6.3E-8 

(Core  mass)/ih  (s) 
c 

2.4E-5 

1.8E-5 

The  first  observation  from  these  results  is  that  the  kinetic  energy 
^p-u>  J /kg  is  small  compared  to  hj^  =  1,43  E8  J/kg,  and  so  was 

indeed  ignorable.  The  next  observation  is  the  tremendous  acceleration 
of  the  propellant  gas  while  being  heated,  although  it  does  not  reach  the 
sonic  speed,  which  is  ah  =  9.  5E5  cm/s,  but  stays  at  very  low  Mach 
n.umber.  Finally,  we  see  that  the  detrained  core  ma.ss  flow,  by  our 
estimate,  could  be  comparable  to  the  propellant  mass  flow.  Further,  it 
empties  the  core  in  a  time  of  order  10"^  seconds.  While  this  estimate  of 
detrainment  is  quite  crude,  and  may  be  on  the  high  side  since  the  diffusion 
coefficient  was  evaluated  at  Th,  still  it  shows  that  detrainment  can  be  a 
big  potential  problem  for  the  stationary-core  concept. 

This  example  calculation  is  meant  to  illustrate  the  sort  of  analysis 
necessary  to  define  the  properties  of  a  stationary-core  laser -heated 
rocket.  Most  of  the  uncertainty  in  the  results  comes  from  lack  of 
detailed  knowledge  of  the  absorption  processes  in  the  propellant,  and  of 


-117- 


the  core  detrainment  process.  Both  of  these  areas  need  further  study 
if  the  stationary -CO re  concept  is  to  be  pursued.  With  some  more  detailed 
models  of  these  two  processes  available,  parametric  studies  could  be 
done  at  various  values  of  Tq,  p^,,  P.  There  may  well  be  preferred 
regions  in  which  to  operate  such  an  engine,  and  parametric  studies 
might  uncover  them. 
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11.  CONCLUSIONS 


Based  on  the  studies  of  laser-heated  rockets  described  in  this 
report,  as  well  as  the  knowledge  gained  in  the  previous  study  reported 
in  Ref.  1,  we  can  draw  some  conclusions  about  the  feasibility  and  opera¬ 
tion  of  these  devices.  We  can  also  identify  important  areas  where  our 
knowledge  is  incomplete,  and  further  research  is  needed  to  put  the  con¬ 
clusions  on  a  sound  scientific  basis, 

1.  The  flowing  core  rocket  starts  with  a  laser  supported 
combustion  (LSC)  wave.  For  a  given  laser  power,  the  cross-sectional 
area  of  the  rocket  engine  is  inversely  proportional  to  the  laser  intensity 
at  which  this  wave  operates.  The  temperature  level  of  the  plasma  core 
and  the  mass  flow  rate  through  the  engine  are  also  determined  by  the 
LSC  wave,  so  a  knowledge  of  its  properties  are  vital  to  the  definition  of 
the  rocket.  However,  there  is  at  present  no  theoretical  or  experimental 
information  on  such  waves  in  hydrogen,  and  even  in  air  they  are  not 
well -under stood.  The  properties  used  here  were  estimates  based  on 
extrapolating  a  simple  air  theory  to  hydrogen.  However,  LSC  wave 
properties  in  hydrogen  represent  the  major  uncertainty  in  the  present 
study,  and  should  be  studied  both  experimentally  and  theoretically  to 
increase  confidence  in  the  feasibility  of  laser -supported  rockets  using 
hydrogen. 

2.  Performance  of  the  flowing -core  streamtube  as  a  rocket  is 
good.  Values  of  Igp  for  expansion  to  1  atm  exhaust  are  1500  to  2500  s 
while  expansion  to  a  vacuum  leads  to  3000  to  3300  s. 

3.  Only  a  limited  number  of  calculations  were  made,  using 
three  laser  powers  at  each  of  two  pressures.  The  LSC  wave  properties 
and  the  core  velocity  distribution  were  not  varied.  To  map  out  the  region 
of  usefulness  of  laser-heated  rockets,  more  extensive  parametric  calcula¬ 
tions  are  needed. 

4.  Breakdown  of  hydrogen  to  produce  an  electron  density 
which  is  comparable  to  that  of  thermal  ionization  can  be  produced  by 
either  an  electric  arc  discharge  or  by  a  focused  laser  beam.  Whether 
this  will  actually  initiate  an  LSC  wave  has  not  been  determined,  and 
probably  needs  to  be  verified  by  experimental  investigation. 
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5.  The  radiation  emission  of  hydrogen  has  been  put  into  a 
form  useful  for  fluid  dynamic  calculations.  The  black  body  portion  of 
the  emission  depends  crucially  on  the  radial  temperature  profiles,  and 
has  not  been  included  in  the  present  one -dimensional  calculations.  Its 
inclusion  would  require  some  model  for  the  radial -variation  of  temperature. 
It  could  be  a  significant  effect,  however,  and  should  be  investigated 
further. 

6.  Based  on  the  transparent  radiation  emission,  radiated  energy 
loss  is  not  a  substantial  fraction  of  the  laser  energy  input,  varying  from 
Z0%  for  smaller  engines  to  4%  for  larger  ones.  This  conclusion  is 
sensitive  to  the  operating  temperature  of  the  plasma  core,  and  could 
change  if  this  temperature  is  even  1000  K  higher  than  calculated  here. 

7.  The  energy  flux  radiated  is  sufficiently  high  to  require  some 
kind  of  active  heat  protection  for  the  walls,  ranging  up  to  75  kW/cm^. 

8.  The  convective  energy  losses  calculated  by  assuming  a 
hydrogen  boundary  layer  at  the  edge  of  the  plasma  core,  without  coupling 
the  loss  to  the  core  flow,  are  too  high  to  be  realistic  for  all  but  the  5000  kW 
laser  power.  A  coupled  calculation,  not  undertaken  here,  is  needed  for 
smaller  power  levels. 

9.  For  the  5000  kW  power,  the  convective  heat  losses  are  of  the 
order  of  30  kW/cm^  near  the  throat,  at  3  atm,  and  perhaps  twice  that  for 
10  atm,  so  they  add  significantly  to  the  need  for  heat  protection,  though 
they  only  occur  in  a  small  region. 

10.  Normal  injection  of  hydrogen  into  the  core  streamtube  seems 
to  be  a  useful  way  of  reducing  the  convective  heat  losses,  based  on  a  few 
sample  csLlculations.  Further  calculations  are  needed  to  establish  this  as 

a  practical  method  of  heat  protection.  Tangential  injection  was  not  studied, 
but  should  also  be  investigated  for  heat  protection. 

11.  First  estimates  indicate  that  particle  absorption  does  not 
seem  a  useful  way  to  protect  the  walls  from  radiation,  because  of  the  high 
mass  loading  of  particles  required.  This  subject  needs  further  study. 

12.  The  length  of  the  plasma  core,  up  to  the  throat,  is  controlled 
by  the  inverse  Bremsstrahlung  absorption  coefficient.  Larger  absorption, 
caused  by  larger  pressures,  leads  to  shorter  lengths.  At  3  atm,  this 
length  is  4  cm,  while  at  10  atm  it  is  0.  5  cm.  These  short  lengths  at  the 
high  pressure  lead  to  short,  stubby  streamtube  shapes,  which  may  be 
undesirable  from  a  fluid-dynamical  point  of  view,  especially  for  larger 
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powers.  This  imposes  a  limitation  on  how  high  a  pressure  level  can 
be  used  for  operation. 

13.  The  stationary- core  rocket  engine  concept  was  studied 
qualitatively  and  semi -quantitatively.  The  physical  processes  involved 
in  this  concept  are  more  complex  than  those  of  the  flowing -core  concept. 
Two  major  uncertainties  were  uncovered.  One  is  a  lack  of  knowledge 
of  absorption  lengths  for  high  temperature  hydrogen  radiation  in  low 
temperature  hydrogen.  The  other  was  a  lack  of  knowledge  of  how  the 
flowing  propellant  gas  entrai  ns  core  gas.  Further  study  of  these  two 
areas  would  enable  estimates  of  stationary -core  rocket  performance  to 
be  made  with  more  confidence. 
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